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abstract 


This report presents an overview of research work focussed on the design and insertion 
of classical models of human pilot dynamics within the flight control loops of V/STOL 
aircraft. The pilots have been designed and configured for use in integrated control system 
research and design. The models of human behavior that llave been considered are: 1) 

A Ki * 

McRuer-Krendel - a single variable transfer function model, 2) Optimal Control Model - a 
multi-variable approach based on optimal control and stochastic estimation theory. These 
models attempt to predict human control response characteristics when confronted with 
compensatory tracking and state regulation tasks. 

An overview, mathematical description, and discussion of predictive limitations of the 
pilot models is presented. Design strategies and closed loop insertion configurations are 
introduced and considered for various flight control scenarios. Models of aircraft dynamics 
(both transfer function and state space based) are developed and discussed for their use in 
pilot design and application. Pilot design and insertion are illustrated for various flight 
control objectives. Results of pilot insertion within the control loops of two V/STOL 
research aircraft ( IjSikorski Black Hawk UH-60A, 25^McDonnell Douglas Harrier II AV- 
8B') are presented and compared against actual pilot flight data. Conclusions are reached on 
the ability of the pilot models the adequately predict human behavior when confronted with 


similar control objectives. 
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I. INTRODUCTION 


The practical problems associated with aircraft flight control system design and 
evaluation are complex and wide ranging. During the design phase, control system 
designers arc faced with the selection of the control parameters that will best fit the system 
performance and control objectives, the vehicle configuration, and the particular situation. 
Powered-lift V/STOL aircraft pose unique problems in controls design since the propulsion 
controls are an integral part of the overall flight control system. Control system response 
characteristics and thus the control parameters are typically chosen from analytic 
evaluations and iterative design methods. Mathematical models of the vehicle and control 
system are developed. Control parameters arc selected by various techniques to achieve the 
desired response characteristics. Evaluations of the control system performance are carried 
out in computer simulation environments or by a closed form approaches. Test input 
deflections (i.e. steps, ramps, disturbances, etc.) are injected into the various control 
mechanisms to directly excite specific closed loop dynamics. Measurement and subsequent 
evaluation of the closed loop vehicle and control system responses determine if the desired 
characteristics have been achieved. 

The analysis of the control system and aircraft dynamics, in this typical design 
approach, provides the designer with valuable insight to the system's general/functional 
operation, but with only a limited basis to gauge the final selection of the control 
parameters. These limitations are associated with the lack of the total system response (i.e. 
the summation of the aircraft, control system, and the human pilot's dynamics). When an 
actively participating human is introduced to the flight control environment, he brings with 
him a complex array of control responses, that serve as inputs to the vehicle and control 
system. These inputs differ significantly from the test inputs used during the initial design 


1 



process. Thus the pilot's dynamic behavior must be carefully considered throughout the 
design and evaluation process. 


When actively involved and participating in a flight control environment, the pilot acts 
as a controls integrator by performing a variety of control tasks via the manipulation of the 
multiple cockpit control mechanisms (e.g. longitudinal and lateral stick, rudder pedals, 
throttle, etc.). In addition, some flight scenarios require the pilot to perform tasks other 
than the control commitments (e.g. communications, navigation, weaponry system, etc.). 
Because the pilot is such an integral component of the overall flight control structure, pilot 
opinion tends to be the major criterion for deciding whether control system performance is 
satisfactory. 

Current approaches to incorporating human response characteristics during the control 
system design involve the use of fixed-base piloted simulators. This phase is a necessary 
step in the overall system design process and provides the ultimate source of human 
response characteristics for evaluation purposes. In addition, the pilot s opinion can be 
directly incorporated as the final gauge of the overall system performance. The direct 
incorporation of fixed-base human piloted simulation within the design phases provides a 
safe and effective environment for controls design, but suffers from problems of cost, 
scheduling, and inconvenience. These stem from the necessity that the development take 
place at remote sites within the simulator structure and schedule. To further complicate 
matters, one must also take into account the acquisition of qualified pilots. 


An alternative approach, to the direct incorporation of actual human pilots, is the use of 
analytic models of human behavior within a computer simulation environment. These 
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"paper pilots" attempt to simulate various aspects of a human's dynamic response 
characteristics when confronted with certain classes of closed loop control objectives. 
Analytic pilot models provide the ability to analyze and evaluate the fully integrated, total 
system response characteristics, before the control system is introduced to a remote manned 
simulator facility. This approach represents a significant advantage in cost, time, and 
convenience by allowing the base line control system design to be completed at the home 
institution and then to be thoroughly tested and adjusted at a remote, manned simulator site. 

Human pilot behavior is, however, very complex. Analytic models tend to be limited in 
their abilities to precisely predictor human behavior in a given situation. These limitations 
stem from the model's inability to fully simulate the human's methods of deriving 
information from a variety of sources (e.g. visual, auditory, etc.), the human's complex 
information processing activities, and his physical methods of applying the control 
commands. In addition, individual human pilots perform/act differently when confronted 
with identical control objectives. Thus, human pilot models can only attempt to represent 
human behavior is a very general sense. 

I.A. A Piloted Control Simulation Structure 

Pilot activities within a flight control environment are directed at a wide variety of 
operations and objectives. The pilot must provide the necessary control, stabilization, 
guidance, and navigation, along with any additional tasks associated with a specific 
mission. The pilot supplies the controls needed to achieve the mission objective by actively 
analyzing his environment and instituting the appropriate control commands. The pilot 
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therefore functions as a control integrator, by acting as the center piece of the entire control 
structure and actively participating in the closed loop control efforts. 

A human pilot flies an airplane by a feedback method. He senses by sight or feels by 
"the seat of the pants", the motion of the aircraft, and manipulates various cockpit control 
mechanisms to minimize the error between the actual and some desired motion. In other 
words, the motion of the aircraft is perceived, both directly and indirectly through the 
visual inspection of cockpit flight instrumentation (e.g. altimeter, artificial horizon, etc), 
external visual cues, auditory (the manner by which a car driver can shift gears by the 
sound of the engine), and by physical means. Through his computational mechanisms and 
thought processes, the pilot assesses the perceived vehicle attitudes/orientations, and 
determines the necessary compensative corrections. The pilot applies the corrections to the 
cockpit control mechanisms, through the physical movement of his muscular system. This 
is a form of negative feedback control, where the controller (pilot), must close the loop to 
achieve some desired, overall control objective. 

To simulate active pilot participation within a closed loop control environment, the 
control structure illustrated in Figure I.A.-l. has been considered. Within this structure, the 
pilot appears as a cascaded compensator that is driven by command based vehicle attitudes 
and orientations from some type of guidance or navigation process. The pilot attempts to 
orient the vehicle' in the manner specified by the command via the manipulation of the 
appropriate cockpit control mechanisms. The control configuration shown in Figure I.A.-l 
assumes that the pilot feedback is based on visual assessments of cockpit instrumentation 
and external visual cues. The assumption of visual feedback reduces the effectiveness of 
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the simulated pilot because of the inherent limitations due to the lack of other forms of 
feedback (i.e. auditory, physical, etc.). 

The research work that is presented in this report is focussed on the development of 
analytic human pilot models, "Paper Pilots", to serve as design tools for controls systems 
research. The models are tailored for use in computer simulation environments involving 
V/STOL research aircraft. The control structure of Figure I. A.- 1 is used as the basis for the 
design, analysis, and insertion of the pilot model within the control loops of the V/STOL 
aircraft that are considered. This report will present an overview of the human models, 
simulated aircraft environments, and inserted pilot results of the research that has been 
conducted. Chapter II presents a description of the models of human dynamic behavior that 
have been utilized. The inherent limitations in their abilities of predicting human behavior 
are discussed. Chapter III discusses the techniques utilized in the design and insertion of 
the pilot models. Aircraft simulation environments and some simplified models of vehicle 
dynamics are introduced for use in the design of the pilot models. Pilot design is 
demonstrated and control objective considerations are discussed. Finally, configurations 
for pilot insertion to the control environment are discussed. Chapter IV introduces the 
concepts involved in the development of the vehicle dynamics models that have been used 
in the design of the pilots. Model structure and parameter identification are discussed and 
considered for their compatibility with the specific pilot configuration. Models of specific 
V/STOL vehicles and regions of the flight envelope are presented. Chapter V illustrates the 
pilot design and control loop insertion strategies for the specific V/STOL vehicles. Pilot 
parameters are presented for the various vehicles and flight scenarios. Results of pilot 
insertion are presented and discussed. Comparisons are made between the pilot models and 
the flight data of actual human pilots performing similar flight control operations. Chapter 
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VI presents a discussion of the results and an overall conclusion of the research that has 
been conducted and poses questions to be considered for future research efforts. Appendix 
A discusses the general mathematical characteristics of the Optimal Control Model of 
human behavior. The continuous time OCM is discussed and a discrete time representation 
is derived. Appendix B presents the derivation of a high order state space model of the 
V/STOL research aircraft (Harrier II AV-8B) for use with the OCM. Appendix C is a users 
guide for the OCM software. The general structure of the OCM code is introduced and 
procedures for installation and application of the OCM are presented and discussed. 
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II. MODELS OF HUMAN PILOT BEHAVIOR 


The human pilot models that have been considered in this research are based on quasi- 
linear models of human behavior in closed loop compensatory tracking and state regulation 
tasks. The predicted/simulated human compensative control responses are generated from 
visual assessments of some displayed error or external visual cue. The models do not 
consider the other techniques that actual humans utilize to obtain information (i.e. auditory, 
physical (seat of the pants)). A human's control characteristics can be simulated by a 
cascade of three linear operators[ 1,2, 3,4,5, 6,7], as shown in Figure II.- 1, and enumerated 
below: 

1. Neuro-Muscular/Motor Dynamics - This operator describes the lags/bandwidth 
constraints imposed by the human's muscular system and is approximated by an 
adjustable, linear, first-order lag given by: 

H»< s > * TJTT <II " 1) 

where T n is the time constant of the neuro-motor response. It is important to note that the 
human's muscular bandwidth is often restricted by the rate limitations of the cockpit control 
mechanism. The time constant of this lag can be selected to accommodate these effects on 
the bandwidth at which the human can exert control operations. 

2. Pure Time Delays - These operators describe various internal time delays associated with 
visual information processing and neuro-motor signal pathways. 

3. Equalization Network - This operator describes the control strategy implemented by the 
human to close the loop in a manner that best fits a given situation. Typically, the human 
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will select the equalization network to provide the dominant closed loop control response 
with a damping ratio (^) in the range, 0.5-0. 8, and a natural frequency (co n ) in the range, 3- 

4 rad/sec. 


The inherent unpredictability of the human's response is simulated by a random 
component called the controller remnant. 

In this research, two types of human dynamics models have been considered: 1) 
McRuer-Krendel - a single variable linear transfer function description, 2) An Optimal 
Control Model - a multi-variable state space approach based on optimal control and 
stochastic estimation theory. Each model is based on differing implementations of the three 
cascaded operator description. This section presents a general overview of the two models, 
their mathematical foundations and their inherent limitations of simulating/predicting human 
dynamic behavior. 

II.A. - McRuer-Krendel Model of Human Dynamics 


The McRuer-Krendel model (MKM) is a single-degree of freedom quasi-linear model 
based on a best fit analysis of experimental human response data [1,2, 3, 6,7]. The model 
uses assessments of visually-based information to produce compensative control 
mechanism displacements. The general form is given by: 


H p (s)=K p e T D s 


(T l s+1) J (T k s+1) 1 

(Tis+i) |(T k s+1) [(s/con^+^n/ojnIs+i] 


(II.A.-l) 
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where H p (s) is the transfer function of the human response, often referred to as the 
describing function, s is the complex Laplace transform variable, the input is the visually 
based error signal, while the output is the corresponding control displacement. McRuer and 
Krendel discuss typical values of the precision model in [1,2,3]. 

Within this model, human activities are represented by physiological and equalization 
sections. Physiologic attributes simulate the limitations and abilities of the human's 
physical mechanisms. The inherent lags associated with the human's visual, information 
processing, and signal transmission systems are modeled by a pure delay. The restrictions 
associated with the muscular system are represented by the system within the brackets. The 
equalization attributes simulate the control strategies employed by the human to achieve the 
required closed loop responses in the form of a lead-lag network. The primary focus of the 
research using this type of human pilot model has been directed toward the fundamental 
control activities of the pilot during various flight control maneuvers. A zero remnant is 
therefore used during the development of the McRuer/Krendel pilot control strategies. 
Remnant selection for these types of piloted flight configurations is described in [8]. 


The model of EQ(II.A.-l) can be simplified to obtain the transfer function, 


Hp(s) 


Krie- T P»(T LS +1) 

(T N s+l)(TiS+l) 


(II.A.-2) 


where very low and very high frequency accuracy is not necessary. This transfer function 
model is illustrated in the simplified block diagram of Figure n.A.-l. The rejection of the 
very low and high frequency content is a reasonable assumption for the human pilot since, 
as discussed before, the bandwidth of the closed loop is 3-4 rad/sec (or 0.48-0.64 Hz). 
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The pure transmission time delay parameter T D is estimated to be between 0.13-0.23 

seconds [1,2, 3, 6, 7]. Although the changes in the time delay can be significant depending 
on the particular control task, the parameter selected at Tq= 0.20 was considered reasonable 
for the purposes of this research. The term l/(T N s+l) is an approximation of the neuro- 
muscular lag of the arm meaning that the pilot can not move his arm faster than the rate of 
this pole. The value of T N is assumed to be constant and approximately 0. 10 seconds. The 
remaining term, Kp(T^s+ 1 )/(TjS+ 1 ), is the equalization part of the model (a time dependent 
variable gain and a lead-lag compensator) whose parameters are altered by the pilot to the 
particular flight configuration and control objective. The constraints on the model 
parameters are as follows: 


0.0£T,<2.50 (T l *T n ) (II.A.-3a) 

0.0<Ti<20.0 (II. A. -3b) 

T n =0.10 OI.A-30 

T d = 0.20 (II. A. -3d) 


The lead-lag compensator/equalizer is based on the assumption that the human is 
required to furnish at least one differentiation and one integration to obtain the desired 
performance, and the constraints on the parameters, T^ and T[ determine how efficiendy the 

integration and differentiation processes are performed by the human. Even though there 
are only a few parameters to be adjusted, the analysis is still not trivial because of the time 
delay, time-varying pilot parameters and time-varying aircraft dynamics. 


II.B. - The Optimal Control Model 


The Optimal Control Model (OCM) [4,5,6, 7] is a multi-variable approach to the 
modeling of human behavior. The OCM is based on the use of modern control and 
estimation theory to simulate/predict human behavior in closed loop control operations. The 
OCM is capable of treating multi- variable systems by incorporating a single conceptual 
framework based on state space techniques. The primary assumption involved in the OCM 
is that a well trained human pilot/operator behaves in a near optimal manner subject to his 
inherent limitations, constraints, and control tasks [4], This optimal behavior is simulated 
by an analogous optimal control system. The optimal control system operates to minimize a 
quadratic performance index in the presence of various system inputs, noises, and 
disturbances. 


The system under control consists of the control element and displays which are 
modeled by a linear state equation and output equation. 


x(t) = A x(t) + B u(t) + w(t) 


(H.B.- 1) 


y(t) = C x(t) 

The "n" dimensional state vector is defined by: 

x = [x 1 ,x 2 ,...,x n ] T 
The human manipulates "m" controls: 

u = [u 1 ,u 2 ,...,u m ] T 


(II.B. -2) 


(II.B. -3) 


(II.B. -4) 


1 1 


and observes "1" system displays (output variables): 


y = [yi»y 2 »— yil T (ii.b.-5) 

It is assumed, from remnant and psychophysical studies of human perception [9] that the 
human can extract position and rate from a single display or external visual cue, but can not 
extract higher derivatives. Thus the output "y(t)" contains those quantities explicitly 
displayed as well as those implicitly derived by the human. This is an important concept 
because it will be directly incorporated in the organization of the vehicle model and the 
strategies associated with displayed information. The disturbance, w(t), is a vector of zero 
mean, white gaussian noise processes and is generally associated with atmospheric 
turbulence when considering aircraft applications. 

The OCM models human behavior in two categories: 1) intrinsic human limitations, 2) 
human control/equalization efforts. Simulation of human limitations is provided by a time 
delay, a neuro-muscular dynamics model, and a controller remnant. The time delays 
associated with visual information processing, neuro-muscular signal propagation, and 
other operations are combined into a lumped equivalent perceptual time delay, T D . It is 

assumed that all outputs are delayed by the same amount. Typically, this delay is on the 
order of 0.2 seconds [4,5,6,7,10,11,12]. Neuro-muscular/motor dynamics are represented 
by an adjustable lag matrix, T N . This lag is not directly modeled as an inherent limitation, 

but is indirectly incorporated by weighting the control rate terms in the cost function of the 
optimal control strategy. The inclusion of a control rate term results in a first-order lag 
being introduced in the optimal controller. This term is utilized to indirectly model the 
physiological limitations of the rate at which a human can perform a control action due to 
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the neuro-muscular/motor dynamics. Controller remnant is modeled by an observation 
noise vector, Vyy(t), and a motor noise vector, V uu (t), where 

E{V yy (t),V w T (o)} = R y 5(t-o) (II.B.-6) 

E{V w (t),V uu T (a)} = Q u 5(t-a) (II.B.-7) 

The observation noise models the inherent uncertainty of the human's visual assessments 
of the displayed information. A separate noise source is associated with each displayed 
output. The noise processes are modeled as an independent, zero mean, white, gaussian 
noise sources. The spectral density is proportional to the mean-squared value of the 
displayed variable, which is basically a signal-to-noise ratio that is on the order of -20dB 

[5.6.7.8.10] , The motor noise models the inherent uncertainty of the human's control 
execution. Like the observation noise, the motor noise is assumed to be independent, zero- 
mean, white, and gaussian. The spectral density is proportional to the mean-squared 
operator output. The motor signal-to-noise ratio is typically chosen near -25dB 

[5.6.7.8.10] . 


The human's equalization network describes the manner in which the human attempts 
to optimize his control strategy to match a given situation. As shown in Figure II.B.-l, the 
human perceives a delayed noisy replica of the system output, y p (t), where: 

y p (t) = y(t-d) + Vyy(t-d) (II.B.-8) 

A 

Estimation of the delayed state vector, X(t-T D ), is accomplished via a Kalman filter. The 
Kalman filter models the human's deduction of the system states from the displayed 

A 

information. A least mean-squared predictor generates a present time state estimate, X(t), 
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A 

from the delayed estimate, X(t-T D ). The predictor models the human's compensation for 
his inherent time delay. The optimal gain matrix, K*, is generated by the solution, in steady 
state, of the optimal regulator problem [13] for the cost function of the form: 



T 



x(t) T Qx(t) + u(t) T Ru(t) + u(t) T Su(t)]dt 

y P (<*) 


} 

(II.B.-9) 


where 

a< t 

Q and R are positive semi-definite 

S is positive definite 

The application of the optimal control, Kalman filter, and predictor require the use of an 
internal reference model (i.e. the model of the vehicle as perceived by the operator) to 
generate their appropriate gains and parameters. Thus the model of the system under 
control plays an important role in the response actions of the OCM. Appendix A.l presents 
a mathematical overview of the OCM in a continuous time representation. For use in the 
simulation model program, the continuous time model is converted to a discrete time 
representation. The concepts involved and the resulting discrete time OCM are presented in 
Appendix A.2. 
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HI. STRATEGIES OF PILOT MODEL APPLICATION 


The inherently different structures of the McRuer-Krendel model (MKM) and the 
Optimal Control Model (OCM) require different approaches in parameter selection and 
control loop insertion. The MKMs arc specifically designed for each control objective and 
region of the flight envelope. The equalization parameters are selected off-line, by Root- 
Locus techniques and arranged in a gain table format. The multi-variable aspects of the 
OCM are directly incorporated in high order control configurations. Cost function 
weightings are selected according to parameters extracted from actual piloted flight data and 
arranged for the specific control situation. This section will discuss the techniques utilized 
to select the pilot's control parameters and the configurations used to insert the pilots within 
the simulated control environment. 

III.A. Aircraft Simulation Environments 


Before introducing the techniques utilized to select and insert the pilot models, it is 
important to discuss the computer-based aircraft simulation environments that were used in 
this research. The principle tools utilized to examine the V/STOL vehicles and evaluate the 
performance of the pilot models were two computer simulation programs [14,15,16,17] 
provided by NASA Lewis. These programs implement nonlinear, total force, large angle 
representations, in six degrees of freedom. These programs provide full flight envelope 
operation and incorporate all on-board stability augmentations systems. 
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m.B. McRuer-Krendel Pilot Models with Static Parameters 


The static MKM pilots are based on the selection of the equalization parameters, for 
specific regions of the flight envelope, through the use of off-line application of Root- 
Locus techniques. In this approach, low order transfer function models of the aircraft 
dynamics are developed for regions of the flight envelope that are of interest. The regions 
of the flight envelope are designated by the vehicle's forward velocity. Pilot equalization 
parameters are selected from the use of a delayed Root-Locus and the direct incorporation 
of the pilot's physically limiting factors. Pilot parameters are then arranged in gain tables 
according to the flight envelope region. This design process is earned out for all cockpit 
control mechanisms (e.g. Lateral and Longitudinal Stick, Rudder Pedals, etc.) and for the 
vehicle attitudes and orientations that are relevant to the flight control objectives. This 
creates a set of Single Variable Pilot Mechanisms (SVPM). The insertion of a static MKM 
pilot is achieved by selecting an appropriate group of SVPMs for the objectives associated 
with the flight scenario. The insertion group has a SVPM defined for each cockpit control 
mechanism and thus has a unique feedback structure. The insertion group is therefore valid 
for only a limited number of flight control objectives. 

III.B.l. - Structure of Aircraft Dynamics Models 


The models of aircraft dynamics for use in the design of the MKM pilots, are based on 
single variable low order, linear transfer functions. The transfer function models are 
derived from the dominant response characteristics of the vehicle dynamics due to the 
injection of test inputs (e.g. impulses, steps, etc.) to the cockpit controls. The dominant 
responses refer to the most pronounced (primary) vehicle attitude reaction due to the 
deflection of a single cockpit control mechanism. The remaining vehicle reactions (i.e. the 
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coupled responses occurring in the other attitudes due to the operation of the cockpit control 
mechanism), are considered secondary. This definition of primary and secondary vehicle 
responses is illustrated in Figure III.B.l.-l. During SVPM development, the secondary 
reactions are ignored and considered as disturbances. The secondary reactions are, 
however, considered in the selection of an insertion group. 

The dominant aircraft dynamics responses were identified by the direct evaluation of the 
time and frequency domain responses of the simulation model programs due to the test 
input deflections of the cockpit control mechanisms. Time based responses were utilized to 
match simplified low order responses. The frequency responses were obtained by Fast 
Fourier Transforms of the time responses. Attempts at the determination of ultra-low 
frequency response characteristics were hampered by the presence of parasitic low 
frequency response modes (Phugoid Modes) and distortions associated with the cross- 
couplings of the secondary variables. Solutions to these problems are considered in 
[19,20]. 

Along with the initial decoupling of the primary and secondary responses, the control 
and response characteristics of each vehicle were separated into two groups: 1) 
Longitudinal Control Set, and 2) Lateral -Directional Control Set. These sets consider the 
effects of the control mechanisms on the overall orientations of the vehicle to the primary 
orthogonal planes [18,19,20]. These sets simplify the selection of the SVPM when 
constructing an insertion group. 
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III.B.2. - Development of SVPM Equalization Parameters 


The development of the Single Variable Pilot Mechanisms are based on the single loop 
control structure shown in Figure III.B.2.- 1. This single variable control configuration is 
organized in such a way that the manipulation of a specific cockpit control mechanism is 
based on the assessment of visual feedback obtained from the observation of a single 
cockpit instrument or external visual cue. The relationship of aircraft dynamics to the 
specific cockpit control is obtained from the analysis of the aircraft dynamics discussed 
previously, and the desired control variable. 

The MKM equalization parameters are selected via Root-Locus techniques. A problem 
that complicates this approach is the pure time delay associated with the human s visual, 
computational, and signal conduction delay model. In the continuous time domain, a pure 
time delay corresponds to an infinite number of poles positioned at S = -°° on the real axis. 
These poles introduce an infinite number of asymptotes that are parallel to the real axis. The 
presence of these asymptotes (specifically the primary asymptote) create significantly 
destabilizing distortions of the asymptotic behavior of the Root-Locus, as shown in Figure 
III.B.2.-2. An additional problem that is presented by the delayed Root-Locus is that many 
of the non-delayed assumptions are no longer valid. 


To overcome the complications associated with the delayed Root-Locus, an 
approximation to the infinite pole set was applied in a non-delayed format. The pure delay 
approximation is given by: 
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(III. B. 2.-1) 


c t d s ~ 

(1 + ^ 

where Tp = 0.2 seconds for this application. 

This approximation not only simulates asymptotic behavior, but also permits all non 
delayed assumptions. To remain within computational limits (Quad Precision), 20 poles 
were placed at S = -100 on the real axis, as shown in Figure III.B.2.-3. The primary 
asymptote of this pole set has a 9° angle of incidence and has an imaginary axis intercept at 
15.84 rads/sec (15.71 rads/sec for the ideal pure delay). This pole set serves as a 
reasonable approximation to the pure delay when considering the closed loop natural 
frequency of the piloted control (W N - 3.5 rads/sec). 


Incorporating the delay approximation and the muscular system limitation pole, a 
modified transfer function for design purposes only) can be given by: 




(10Q) 20 

(S + 10)(S + 100) 



(III.B.2.-2) 


where a £ 0.8 

b > 0.0 

III.B.3. Multi-variable McRuer-Krendel Pilot Insertion Techniques 

When a participating human pilot is introduced into the control loops, he uses all 
cockpit control mechanisms to provide an operational control by employing visual, audio, 
and other forms of feedback cues. For the purpose of simulation, a comprehensive human 
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response pilot model is constructed by integrating various single variable pilot mechanisms 
(SVPM) into a multi-variable structure (an insertion group), as shown in Figure III.B.3.-1. 

The multi-variable pilot structural configuration is defined by the manner in which 
visual feedback is interpreted and applied to the cockpit control mechanisms via a specific 
set of Single Variable Pilot Mechanisms (SVPM). The intrinsic limitations of each 
configuration (insertion set) make it applicable to only a specific set of flight control 
maneuvers. This is primarily due to the limited number of feedback paths that are available 
(i.e. the number of feedback paths equals the number of cockpit control mechanisms). Each 
type of flight maneuver or command sequence is therefore associated with a specific multi- 
variable pilot configuration. In general, a command maneuver will be described by a set of 
vehicle attitudes and/or rates that define the new orientation that the vehicle is required to 
attain. The intricacy of the maneuver defines the number of attitudes and/or rates that are 
simultaneously involved in the operation. 

The flight control objectives associated with simple maneuvers require the control of 
only one primary vehicle attitude or rate. The remaining secondary variables are monitored 
and regulated (regulation set) to preserve the stabilized aspects of the vehicle orientation 
(e.g. level flight). This type of operation can be performed by the configuration shown in 
Figure III.B.3.-2. A regulation set is defined by the single variable pilot mechanisms that 
are associated with the flight control variables needed to maintain a stabilized operation in a 
specific maneuver. The selection of the regulation set is based on the primary command 
attitude and the general operations involved in the execution of the maneuver. Large scale 
secondary responses must be compensated, proper feedback paths must be allocated, and 
the appropriate SVPMs must be utilized. This is similar to the manner in which a well 
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trained pilot reconfigures his feedback for a given situation. An additional concern in the 
selection of an insertion set is the use of multi-function cockpit control mechanisms. These 
controls are utilized differently depending on the control objectives. A good example of a 
multi- function control mechanism is the rudder pedals. During level forward flight, the 
pedals are used to make minor heading corrections, while during coordinated turns, they 
are used to minimize sideslip. This type of control mechanism is handled through the 
assignment of the feedback path and equalization network. 


III.C. - Optimal Control Model 

The multi-variable structure of the OCM makes it's insertion to the control loop 
relatively simple. The application of the OCM is based on parameter selection in two basic 
categories: 1) Control task related, 2) Pilot model related. This section will provide an 
overview of the general construction of the OCM. 

III.C. 1. - Control Task Description 

The first concern is the description of the control task. The control task must be 
described analytically, this includes the specifications of the system under control (vehicle) 
and the objective of the control activity (trajectory). As indicated previously, the vehicle is 
represented by a state space model. Care must be taken in selecting the vehicle model due to 
its implicit use in the formulation of the optimal control gains, Kalman filter and predictor. 
The state and output vectors must be chosen in such a way as to not limit the manner in 
which the OCM will extract estimates of vehicle orientation from the visually displayed 
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information. Control objectives are reflected in the cost function parameters. The 
controller's specific task is to choose a control input, on the basis of observing the 
displays, so as to minimize a weighted sum of the averaged state and control variables. 


Once the vehicle and the control task have been specified, determination of the 
displayed variables is relatively straight forward. The control task, can at times, indicate the 
variables that are considered, or the displayed variables can be concluded from the available 
cockpit instrumentation. In certain control tasks a marked target, in the form of an external 
visual cue, is used. The variables available from the target are therefore related to it's 
markings and are thus described by the display vector associated with the target. As 
mentioned previously, the displayed variables include the quantities explicitly displayed 
plus their first derivatives. 

III.C.2. - Pilot Description 

The OCM pilot is described by four parameters: 1) Time delay, 2) Observation noise, 
3) Motor noise, 4) Neuromuscular lag. The overall structure of the OCM is based on 
optimal control theory, but the theory does not provide the parameter selection. This 
information is typically obtained from human performance data. 

The time delay in simple compensatory tracking and state regulation tasks is generally 
on the order of 0.1 to 0.3 seconds [4,5,6,7,10,11,12]. In complex tasks, the time delay is 
difficult to determine. Values near 0.2 seconds have shown to be reasonable choices from 
human performance data [7,10,12]. 
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Observation noise plays an important role in the estimation problem, because it tends to 
be the dominant source of controller remnant. Various experiments have been performed 
[8] to obtain reasonable estimates of the observation noise spectral density. Typical values, 
for simple tracking tasks are on the order of -20dB while for complex operations -lOdB has 
shown good results [10]. 

Motor noise is a difficult quantity to extract from human performance information. 
Typically, model matching techniques are incorporated to determine reasonable values. 
Spectral density signal-to-noise ratios on the order of -25dB have been indicated for 
relatively simple tasks [5,6,7] and near - lOdB for complex tasks [10]. The effects of motor 
noise, however, do not appear to be great and have even been neglected in some cases [5]. 

The muscular system model of the OCM is based on the subjective weighting of the 
control rate terms of the cost function. The values of the lag matrix, T N , must be chosen for 
the specific muscular activity and cockpit control dynamics. Classical values of T N =0. 1 do 

not appear to be valid when considering complex control tasks or stiff control mechanisms. 
Time constants on the order of T N =1.4 - 0.2 have been shown to more closely agree with 

pilot data [6,7,10,11,12]. In addition, the muscular system involved with the use of the 
legs (for executing control manipulation of the rudder pedals) must take into consideration 
the inherent reduction in bandwidth. 

III.C.3. - Difficulties in applying the OCM 

The OCM provides a great many algorithmic and computational advantages in the 
quantitative estimation of a human operator's dynamics. There are, however, some 
difficulties that arise during the application of the OCM. The first relates to the explicit 
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requirement that the human pilot description be based on an internal model of the human's 
inherent characteristics, the dynamics of the system under control, and external 

A 

disturbances. To provide a present time state estimate, X(t), the system matrices (A,B,C) 
.system disturbances, the human time delay, observation and motor noises must all be 
known. To generate the controller's optimal state variable feedback gains the A and B 
matrices along with the control objective weights of the cost function are required. 
Essentially, this amounts to a complete knowledge of the pilot from man-machine system 
concepts. The OCM requires a very accurate internal model if it is to adequately function in 
a manner consistent with human behavior within the control environment. 

The second difficulty stems from the fundamentally difficult problems associated with 
identifying the pilot's internal model parameters from experimental data. In addition, the 
optimal control strategy suffers from a degree of over parameterization. From an 
identification viewpoint, observation and motor noises are unresolvable and optimal control 
and state estimation gains can only be determined from the matching of experimental data or 
through the similarity transformation of the model [6]. 

The final problem is associated with the specification of the control objective cost 
function. The cost function parameters must be selected in accordance with the control task 
and thus the OCM designer must speculate on the parameters that will be of importance to 
the actual human pilot 
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IV. MODELS OF V/STOL RESEARCH AIRCRAFT 


Two types of V/STOL research aircraft have been used in this study: l) Sikorski Black 
Hawk UH-60A - a modem high performance helicopter and 2) McDonnell Douglas Harrier 
II AV-8B - a thrust vectored jet fighter. Aircraft dynamics where simulated by nonlinear 
computer program models [14,15,16,17], provided by NASA-Lewis. Both programs 
implement total force, large angle, nonlinear representations of the individual aircraft 
dynamics in time based computer simulation environments. These vehicle definitions 
provide full flight envelope operation and support the onboard flight stabilization and 
control systems. 

The vehicle models that have been developed in this study are designed to complement 
the structures of the individual pilot models. The McRuer/Krendel pilots require the use of 
low order/decoupled transfer function vehicle models , while the OCM model relies on high 
order/coupled state space representations. Linearized models of the aircraft dynamics (for 
use in the design of the pilot models) where developed by a mix of analytic models of 
vehicle motion and direct analysis of the time and frequency responses of the simulation 
program to test deflections of the cockpit control mechanisms. The following sections 
provides an overview of the vehicle, the model structures, parameter identification 
techniques, and the resulting vehicle models. 

IV. A. Harrier II AV-8B : A Thrust Vectored Jet Fighter 


The Harrier AV-8B's thrust vectoring capabilities make it a truly unique aircraft. The 
regions of the Harrier's flight envelope can be defined by the direction of the thrust vector, 
and consists of the three basic flight configurations shown in Figure IV. A.- 1. The high 
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speed configuration, shown in Figure IV.A.-l.a, is characterized by the thrust vector being 
directed forward (i.e. nozzle jet vectors aft). In this mode, the propulsion system supplies 
the forward thrust component in a manner that is common to conventional aircraft. The lift 
and control components are supplied by the aerodynamic surfaces as they are forced 
through the atmosphere. The magnitude of the thrust vector supplied by the propulsion 
system is used to control vehicle speed. The stabilator is used to control the angle of attack 
and altitude. 

The transition mode, shown in Figure IV.A.-l.b, is described by a general loss of 
aerodynamic responsiveness. As forward speed decreases, aerodynamic surfaces loose 
their ability to provide necessary lift and control functions. As the name implies, the vehicle 
control actions are in a transition between atmospheric flight and powered-lift activities. In 
general, sustained flight in this region is avoided by the typical maneuvers associated with 
acceleration to and deceleration from the high speed envelope. 

The low speed configuration, shown in Figure IV.A.-l.c, is characterized by the thrust 
vector being directed upward (i.e. nozzle jet vectors downward). The lack of sufficient 
forward velocity requires that the propulsion system provide the lift components (powered- 
lift). A closer examination reveals that the propulsion system supplies the forward thrust 
and the primary lift component in a manner similar to a helicopter's main rotor. The 
magnitude of the thrust vector is primarily utilized to control the altitude and it's direction is 
used to adjust forward speed. 

The interesting region of the Harrier's flight envelope occurs during low speed 
powered-lift activities. In this region, the aerodynamic properties of the aircraft 
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significantly differ from those of high speed conventional flight. During low speed flight 
the components of lift produced by aerodynamic means are small. The vehicle relies 
primarily on lift components supplied by the propulsion system. In addition, the aircraft's 
aerodynamic control surfaces (e.g. rudder, ailerons, etc.) no longer function as the primary 
control mechanisms. The Harrier relies on the Reaction Control System (RCS) to provide 
the additional control components that are needed to maneuver the aircraft. 

The primary focus of the research on pilot models for the Hamer II AV-8B has been 
directed toward the low speed - powered-lift region of the flight envelope. Dynamic 
response tests and general use of the simulation model programs were conducted while in 
trimmed forward flight at speeds ranging from hover to 35 knots. The vehicle was 
configured with the landing gear down, flaps extended to 60 degrees, and the lift 
enhancement devices fully extended. The SAS was enabled to provide damping of angular 
rates. 


IV.A.l. - Harrier Control Structure 

When operating in the low speed region of the flight envelope, the Harrier's control 
structure shows signs of a modest decoupling of the longitudinal and lateral dynamics. 
Longitudinal and lateral control sets are defined by the control strategies that are associated 
with the cockpit control mechanisms. The longitudinal control set specifies the control 
mechanisms and their associated reactions that primarily influence vehicle responses in the 
longitudinal plane (X-Z body plane). The three cockpit control mechanisms that operate 
within this region are: 1) nozzle angle control, 2) throttle, 3) longitudinal stick. 
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Nozzle Angle Control 


The nozzle angle controls the direction of the propulsion system thrust vector. A typical 
low speed - powered-lift thrust diagram of a decrease in nozzle angle is shown in Figure 
IV.A.l.-l. The nozzle angle for this region of the flight envelope is large. The control 
structure is very similar to that of the longitudinal cyclic of a helicopter. Minor changes in 
the direction of the thrust vector will tend to dominate the forward thrust component. This 
indicates that the nozzle angle will dominate the control of the vehicle’s forward speed. 
Small variations of the nozzle angle in low speed flight will have only a small effect on the 
primary lift components. The moment arm associated with the thrust vector displacement 
from the vehicle's center of gravity will induce a relatively small pitching torque in the 
longitudinal plane. 

Throttle 

The throttle controls the magnitude of the propulsion system's thrust vector. A thrust 
diagram of an increase in engine speed is shown in Figure IV. A. 1.-2. This control 
behavior is similar to the main rotor collective of a helicopter. Changes in the magnitude of 
the thrust vector will tend to dominate the lift component and thus the vehicle altitude. 
Minor changes in the thrust vector magnitude will have only a small impact on the velocity 
components and the thrust vector displacement from the vehicle center of gravity will 
induce a relatively small pitching torque in the longitudinal plane. 


Longitudinal Stick 


The longitudinal stick controls the stabilator angle of attack and the vents of the RCS's 
forward and aft jets. During powered-lift activities the auxiliary thrust components and the 
associated moment arms (due to their physical configurations) of the forward and aft RCS 
jets to induce pitching torque responses about the Y body axis in the longitudinal plane. 
Figure IV. A. 1.-3 illustrates the auxiliary thrust components produced by deflecting the 
longitudinal stick. Engine air bleed of the RCS will tend to effect the propulsion system 
performance. The aircraft longitudinal inertial components tend to induce sluggish 
responses. During low speed flight, the longitudinal stick is primarily used to control 
longitudinal orientations (e.g. angle of attack, pitch angle). It is important to note that the 
primary thrust vector of the propulsion system will be redirected during pitching maneuvers 
if the nozzle angle is fixed. This can result in changes in forward speed and altitude. 


The lateral - directional control set specifies the control mechanisms and their associated 
reactions that influence vehicle responses in the lateral plane (Y-Z body plane) and in the 
directional plane (X-Y body plane). The two cockpit control mechanisms that operate in 
this set are: 1) lateral stick, 2) rudder pedals. 

Lateral Stick 

The lateral stick controls the ailerons and the vents of the wing-tip RCS jets. When 
operating in the low speed region of the flight envelope the auxiliary thrust components and 
the moment arms of the wing tip RCS jets are used to induce rolling torques about the X 
body axis. Figure IV. A. 1.-4 illustrates the auxiliary thrust components produced by 
deflection of the lateral stick. Again, the air bleed of the RCS jets will effect engine 
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performance. The lateral inertial components and the RCS wing tip geometry create a 
highly responsive rotational moment about the X body axis. The rolling motion redirects 
the lateral components of the thrust vector and will therefore tend to dominate the lateral 
velocity components. The lateral plane components tend to be coupled to reactions in the 
directional plane. 

Rudder Pedals 

The rudder pedals control the rudder angle and the vents of the tail-end RCS jets. 
During powered-lift operations, the thrust components and moment arm of the tail-end 
RCS jets are used to execute yawing and lateral control maneuvers. Figure IV. A. 1.-5 
illustrates the auxiliary thrust components produced by deflection of the rudder pedals. 
Again, RCS air bleed will effect engine performance. The RCS thrust vectors induce 
rotations about the Z body axis. This allows rudder pedal control of the lateral velocity 
during turn coordination tasks (sideslip reduction) and heading regulation. 

IV.A.2. - Low Order Transfer Function Models of Harrier Dynamics 

This section presents the development of a set of low order linear transfer function 
models (for use in the design of the McRuer/Krendel pilot models) that are based on an 
investigation of the fundamental behavior of the nonlinear aircraft simulation model. The 
low order models describe dominant, decoupled, short period vehicle dynamics that are 
induced by deflections of specific cockpit control mechanisms [20]. These models do not 
attempt to account for the inherent couplings of the vehicle dynamics or any long term 
response characteristics (e.g. longitudinal phugoid modes, lateral spiral or dutch roll 
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modes). Long term dynamic modes and cross couplings are neglected and considered 
secondary to the dominant responses* 

The low- order linear transfer function models of the aircraft dynamics where obtained 
by an analysis of the time and frequency responses of the simulation program to test 
deflections of the cockpit control mechanisms. To simplify and limit this discussion, 
several examples of low order model identification will be presented. A comprehensive 
analysis and evaluation of the transfer function models can be seen in [20]. 

Pitch Angle/Longitudinal Stick Model 

As shown in Figure IV.A.1.-3, the longitudinal stick dominates the pitching motions. 
The time based, short period pitch rate response while in a near hover, due a 1 inch impulse 
deflection of the longitudinal stick can be seen in Plot IV.A.2.-1. This type of time 
response can be modeled by a simple pole residing on the negative real axis of the 
Laplacian complex plane. A low order transfer function model that describes this type of 
time response is given by: 

q(S) Kthetdt (IV.A.2.-1) 

5 e (S) S + a THETA 

where 

9 

K thetdt = 0.71-^_j (IV.A.2 

thetdt inch-sec 2 

a THETA = sec 1 (IV.A.2 

The pitch angle model is the direct integral of the pitch rate model and is given by: 


.-2a) 

■-2b) 
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6(S) K thet (IV. A. 2. -3) 

5 e (S) S(S + a T HETA^ 

The parameters of these equations vary only slightly, due to the relatively small changes in 
aerodynamic effects within the low speed region of the flight envelope. 

Roll Angle/Lateral Stick Model 


The lateral stick's dominant effect on rolling motions can be seen in Figure IV. A. 1.-4 
tends to dominate the rolling motions. The time based response of the roll rate due to a 1 
inch impulse on the lateral stick is shown in Plot IV. A. 2.-2. The roll rate response is 
similar to the pitch rate response, but has a much shorter time constant. The transfer 
function model is given by: 

c KpHTDT (IV.A.2.-4) 

8 a (S) s + a PHi 


where 


K — 1 RQ 

Kp HIDT" L8y inch - SeC 2 
a PHj = 3.45 sec -1 

The roll angle model is the direct integration 


(IV. A. 2. -5a) 
(IV. A. 2. -5b) 

of the roll rate model and is given by: 


«§> fPHTDT (IV.A.2.-6) 

5 a (S) S(S + a PHF 


Parameter variations are insignificant within the low speed environment. 
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Forward Velocity/Nozzle Angle Model 


As a Final example of the low order transfer function models, consider the response 
characteristics of the nozzle angles effect on forward velocity. The time based response of 
the forward velocity due to a 5 degree impulse of the nozzle angle is shown in Plot 
IV. A. 2.-3. The short period response can be modeled as a step function (i.e. integral of the 
impulse input). The ramping response appearing in the later phases of Plot IV. A. 2. -3 is 
due to the long term phugoid effects. The transfer function model of an integrator is a 
simple pole residing at the origin of the complex Laplace plane, and is represented by: 


k vel 

&r(S) 


(IV.A.2.-7) 


where 


k vel ” 


-0.017 


knots 

deg-sec 


(IV.A.2.-8) 


IV.A.3. - High Order State Space Models of Harrier Dynamics 

The high order, coupled state space vehicle models that have been developed for use in 
the design of the OCM pilot's internal reference model, are based on a set of generalized 
linear, first-order differential equations [21,22], that describe the vehicle motion. The 
equations of motion are of the form: 

^(t) = A x(t) + B u(t) (IV.A.3.- 1) 
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The state vector, x(t), represents the perturbations from trim of the vehicle’s pseudo-body 
axis variables. To maximize the overall usefulness of the internal reference model and to 
apply the control tasks (described in later sections), the state vector of Table IV. A. 3.-1 was 
used. Utilizing this large order state vector provides for a flexible internal reference model. 
The use of pseudo-body axis variables results from the manner in which the human pilot 
model will interpret the flight control environment via a mix of external visual cues and 
instrumental feedback. A close examination of the state vector variable selection and 
organization shows that the state vector is made up of a set of vehicle body angles and 
positions, along with their first derivatives (i.e. angular rates and body velocities). 

The control vector, u(t), represents the deviations from the trim positions of the cockpit 
control mechanisms and is defined in Table IV. A. 3. -2. The use of the cockpit control 
mechanisms is due to the manner in which the human will institute his control actions upon 
the vehicle. 

The output or measurement equation is of the form: 

y(t) = C x(t) (IV.A.3.-2) 

This equation provides the relationship of the variables that are displayed to the pilot from a 
linear combination of states. The structure and organization of the state vector permits the 
displays to primarily include vehicle positions and body angles. This is indeed a good 
structure because many of the cockpit displays and information available from external 
visual cues are in the form of a position indication. It is assumed that the pilot will therefore 
derive the first derivative information from the displays, and thus the entire state vector can 
be estimated given the proper display organization. 
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In the most general aircraft model, the elements of the system matrix A, and the control 
distribution matrix B, consist of two basic types. The first consists of inertial and 
gravitational components that are obtained analytically from the equations of motion. The 
second consists of partial derivatives associated with aerodynamic forces and moments. 
Due to the use of the low speed region of the flight envelope, many of the aerodynamic 
terms can be neglected. In addition, no attempt will be made to incorporate the SAS, 
instrumental, actuation, or cockpit control mechanism dynamics as components external to 
the state model. Instead, these dynamics will be incorporated directly within the state 
model. 

The total linearized vehicle dynamics can be described by a completely coupled state 
model given by: 


f— * 

ii 

^long 

A 

^Mat-long 

x(t) + 
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T5 
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A 

_ Iong-lat 
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-®long-lat 
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u(t) 


(IV.A.3.-3) 


In general, the cross coupling terms, A long . lat , A lat . long , B long _ ]at , and B lat . long , can be 
ignored because of their limited secondary response characteristics and the low speed 
assumption [21,22,23]. Extensive testing of the simulation model program showed this to 
be true (see Appendix B). It is interesting to note that the lateral system did not show any 
signs of response excitation due to longitudinal activities, which results in: 

A = b . =0 (IV.A.3.-4) 

^long-lat “ ^long-lat “ 
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It is believed that this response characteristic is due to the absolute symmetry of the 
computer simulation model. The lateral system components did, however, induce relatively 
small reactions within the longitudinal system, thus: 

^lac-long ~ ®lat-long ~ ® (IV. A. 3. -5) 

The effects of the lateral components will be discussed in the development of the 
longitudinal model in Appendix B. 

Longitudinal Dynamics Model 


A generalized linear representation of the core longitudinal dynamics is given by: 


V u 


X u 

X 9 

X q 



V u 


x e 

X 1 

x T 

e 



0 

0 

1 

0 


e 

+ 

0 

0 

0 

q 


Mu 

0 

M, 

Mw 


q 


Me 

Mj 

Mj. 

-v w _ 


-Zu 

Ze 


z w - 


- V w- 


-Z e 

z j 

Zj- 


6) 



(IV.A.3.- 



The core dynamics are those variables which are the principle components involved in the 
description of the motions and responses. The remaining variables are the direct integration 
of the core dynamics. Using various low speed assumptions (see Appendix B), a variety of 
simplifications were made and verified by an analysis of the responses of the simulation 
model program. The longitudinal system of EQ(IV.A.3.-6) was simplified to: 
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(IV.A.3.-7) 

This model describes the generalized low speed dynamics of a thrust vectored aircraft. The 
response characteristics of this state model were examined and compared against the 
classical V/STOL responses [21,22] and those of the Harrier simulation program 
[15,16,17] (see Appendix B). The comparisons showed and explained many interesting 
response modes associated with the Harrier simulation programs operations. An example 
of this, can be seen in the long term pitch rate response due to an impulse on the 
longitudinal stick while operating outside the ground effects region, as shown in Plot 
IV. A. 3.-1. The low frequency oscillatory response characteristics (noticeable in the latter 
phases of the response) can be attributed to the Phugoid mode [21,22] (see appendix B). 
The Phugoid response characteristic rarely troubles pilot activities because of it s ultra-low 
frequency content. This is similar to driving a car that is "out-of-alignment" . Tlte car driver 
simply compensates by providing an offset at the steering wheel that is necessary to 
overcome the misalignment. A closer examination of the system parameters that introduce 
the Phugoid response characteristics revealed that certain parameters could be neglected 
when operations are primarily directed at the pilot frequencies (i.e. pilot operations are 
directed at the short period dynamics). Applying these additional simplifications to the 
model of EQOV.A.3.-7) resulted in the low speed/powered-lift longitudinal model for pilot 

frequencies as given by: 
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This model represents the core longitudinal dynamics that were incorporated in the OCM's 
internal reference model. 


Lateral-Directional Dynamics Model 


A generalized linear representation of the core lateral-directional dynamics is given by: 
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Applying the low speed assumptions of Appendix B, a variety of simplifications were 
obtained, which resulted in the following low speed/powered-lift lateral-directional model. 
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(IV. A. 3. -10) 
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The response characteristics of this model were examined and compared against the non- 
linear simulation program and the classical V/STOL responses (see Appendix B). As in the 
longitudinal case, many of the low frequency lateral modes could be neglected (e.g. spiral 
and dutch role modes). Applying these further simplifications to the model of EQ(IV.A.3.- 
10), a low speed/powered-lift lateral-directional model for pilot frequencies is given by: 
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This model represents the core longitudinal dynamics that were incorporated in the OCM's 
internal reference model. 


IV. B. Black Hawk UH-60A: A High Performance Helicopter 


The primary focus of the research on pilot models for the Black Hawk helicopter was 
directed at developing a group of McRuer-Krendel pilots that spanned the flight envelope. 
This section discusses the development of a set of low order linear transfer function models 
that where obtained by an analysis of the nonlinear simulation model responses due to 
cockpit control mechanism operations. To simplify and limit this discussion, only a few 
transfer function models will be examined. The full set of vehicle models obtained are 
presented in [19]. 
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During the initial stages of the pilot model development a linear point mass, small 
perturbation, state space model [24] of Black Hawk dynamics was utilized. This model 
was linearized about trimmed flight conditions and considered only the pure body 
dynamics. Augmenting this model with actuation and automated systems created a model 
whose complexity approached that of the nonlinear simulation model. The high order of 
this model created significant difficulties in the design and evaluation of the low order pilot 
models. For this reason, low order transfer functions where utilized to create more 
understandable vehicle models. As suggested above, these transfer function models where 
developed and used in a manner similar to the low order Harrier models. 

Simulation model response tests where conducted while in trimmed flight at forward 
velocities of 20, 40, 60, 80, and 100 knots. The initial test scenarios where carried out with 
all onboard automatic control systems disabled and showed low frequency divergence. To 
improve the overall response characteristics, the pitch bias actuator, automatic tail stabilator 
control, and stability augmentation systems (both digital and analog) where enabled. The 
results showed substantial improvement. 

IV.B.l. - Control Structure and Low Order Transfer Function Models 

The cockpit control mechanisms and their associated vehicle responses where 
decoupled into the longitudinal control set (Longitudinal cyclic stick (S lo _ c ), Main rotor 
collective stick (5 mr c ) and the lateral control set (Lateral cyclic stick (8| a ^ c ), Tail rotor 
collective pedals (5^ c )). A more detailed analysis of the control sets and their associated 

response characteristics is presented in [19]. The following model identification examples 
illustrate the fundamental approach utilized. 
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Pitch Angle/Longitudinal Cyclic Stick 


Consider the pitch angle dynamics associated with the manipulation of the longitudinal 
cyclic stick. Plot IV. B. 1-1 and Plot IV.B.1-2 show the pitch rate responses due to an 
impulse on the longitudinal cyclic stick while flying at 60 and 100 knots, respectively. Both 
plots show a damped sinusoidal response characteristic whose damping ratio decreases 
with increasing forward velocity. This type of response can be approximated by the 
transfer function shown below. 


-9 (IV.B.l.-l) 

5 l0 .c S 2 +25 0 w 0 S+w 0 2 


The pitch angle response is given by: 


9 Kj 

Vc S(S 2 +25 0 w 0 S+w 0 2 ) 


(IV. B. 1.-2) 


Where the pitch rate is measured in radians/second and the pitch angle is measured in 
degrees. Evaluating the pitch rate response throughout the flight envelope resulted in the 
parameter values listed in Table IV.B.l.-l. 


Roll Angle/Lateral Cyclic Stick 

Helicopter roll dynamics are dominated by the operation of the lateral cyclic stick. An 
analysis of the response characteristics of the simulation model program indicated slightly 
different roll reactions, one for low speeds (20 to 40 kn) and one for high speeds (60 to 
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100 kn). This type of behavior has been attributed to the low/high speed mode switching of 
the on-board yaw SAS near 60 kn. The low speed model is given by the transfer function: 


(j> 820 

^la_c < ‘ 20 ’ 40 ' ) S(S 2 +25pW p S+Wp 2 ) 2 


(IV. B. 1.-3) 


where 5 p ~ 0. 15 and w p ~ 9.0 rads/sec. 
The high speed model is given by: 


4) __ 210(S + 3) 

5, a c (6a80,100) ~ S(S 2 +25pWpS+w p 2 ) 2 


(IV. B. 1.-4) 


These models do not indicate the resonant behavior that can be seen in the time and 
frequency responses. Figure IV.B.1.-3 shows the frequency response of the roll rate at 80 
kn. A relatively large peak-notch type characteristic can be seen near 2 and 22 rads/sec. The 
high frequency resonance is attributed to the main rotor. The low frequency resonance is 
associated with aerodynamic phenomena. The low frequency peak is, however, 
troublesome because the pilot will attempt to close the control loop near this frequency. 

Altitude/Main Rotor Collective 


The altitude components of the helicopter's flight dynamics are dominated by the main 
rotor collective. Figure IV.B.1.-4 shows the altitude rate response due to the injection of a 
1 inch step on the main rotor collective stick while at a traveling at 60 kn. This type of First 
order response suggests the transfer function model given by: 
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5 mr_c S+aALT 


(IV. B. 1.-5) 


The altitude response is therefore given by: 


ALT __ K-alt 

5mr_c S(S+aALl) 


(IV. B. I .-6) 


Evaluating the altitude rate response throughout the flight envelope resulted in the parameter 
values listed in Table IV. B. 1.-2. 
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V. PILOT DEVELOPMENT AND INSERTION 


This chapter will illustrate the development and insertion of a variety of pilot models 
within the control loop of the non-linear simulation programs. First, a set of McRuer- 
Krendel pilot mechanisms will be developed for the Harrier and the Black Hawk. A group 
of flight control maneuvers is defined and insertion set configurations are organized. 
Results of multi-variable pilot insertions are illustrated by examining several examples of 
Harrier and Black Hawk flight control maneuvers. Finally, the OCM is developed for use 
in the Harrier. Control tasks and display configurations are discussed. Task and pilot 
parameters are chosen. Results of OCM insertion are examined and compared to actual 
human pilot flight data. 


V.A. - Static McRuer-Krendel Pilots 

The concepts involved in the selection and insertion of the McRuer-Krendel model are 
based on the selection of a set of single variable pilot mechanisms (SVPM), to provide a 
variety of control functions associated with each of the cockpit control mechanisms. From 
this set/pool of SVPMs, a group (insertion set) is chosen, one SVPM for each cockpit 
control mechanism. The insertion set is chosen from an analysis of the control functions 
associated with the expected flight maneuvers. This section will illustrate the design of 
some example single variable pilot mechanisms and the selection of various insertion sets 
based on the flight control objectives. Examples of insertion set operation within the 
simulation environment of the non-linear program are presented and discussed. 
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V.A.l. - Design of the Single Variable Pilot Mechanisms 


This section presents some examples of the design of the single variable pilot 
mechanisms (SVPM). The objective is to develop a set of SVPMs, where each pilot 
mechanism is selected to perform a specific control function. The elements of the multi- 
variable configurations will ultimately be chosen from this pool of SVPMs. 

During the design of the individual pilot's, attempts were made to maximize the system 
TYPE. This simulates the pilot's ability to choose the control parameters that will tend to 
improve his error tracking ability. The ability to obtain the desired closed loop response 
characteristics was at times hampered by control mechanism restrictions and difficult 
vehicle response modes (e.g. high frequency complex poles near the right half plane, 
resonant behavior near the closed loop frequencies). Gain limitation and reduced closed 
loop bandwidths were required, in some cases, to keep pilot responses within limits. 

Harrier Pitch Control Pilot 

This pilot provides the position control of the pitch angle by considering visual 
assessment of the pitch angle and manipulating the longitudinal stick. Typically the human 
pilot would derive this feedback from the artificial horizon or from some type of external 
target. This individual SVPM will be used to maintain level flight or to perform specific 
pitch angle positioning. In addition, the pitch angle control can be applied to velocity 
control when performing fixed nozzle angle maneuvers. 

Returning to the pitch angle response model of EQ(IV.A.2.-3), one can see that the 
integrating pole of the pitch angle response provides an intrinsic TYPE I system 
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characteristic. Attempts at introducing the pilot's equalization pole as a second integrating 
pole (to achieve a TYPE II system characteristic) did not provide adequate response 
characteristics due to the positioning restrictions of the equalization zero (a ^ 0.8). The 
pilot's equalization profile was relaxed and the equalization pole was removed to S = 8. 
The equalization zero was placed at S = 3 to maximize the Root-Locus dominant pole 
placement characteristics and provide the necessary phase considerations. This configured 
the equalization as a lead network which in-tum improved the low frequency response 
characteristics and oriented the pilot to be error rate sensitive. The resulting pilot transfer 
function is given by: 


^HARRIER,. 
0 THETA (S) 


4.3 (S + 3) e' 0,2S 
“ (S + 8)(S + 10) 


(V.A.l.-l) 


The Root-Locus for this configuration is shown in Figure V.A.l.-l. The dominant 
poles were placed near the desired locations (w n = 3.1 rads/sec, £ = 0.74). 

Additional Pilot Mechanisms for the Harrier 


Table V.A.l.-l lists the remaining single variable pilot mechanisms for the Harrier, that 
were developed in [20]. The models are based on the form: 


~ HARRIER, x 
0 PILOTS 


K (S + A) e' 0-2S 
~ (S + B)(S + 10) 


(V. A. 1.-2) 
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Black Hawk Pitch Control Pilot 


As in the Harrier pitch control pilot, this pilot also controls the pitch angle position by 
visually assessing the pitch angle and operating the longitudinal cyclic stick. The different 
response characteristics of the Black Hawk required a significantly different pilot 
configuration. The transfer function model of the Black Hawk's pitch dynamics, 
EQ(FV.B.l.-2), permits the construction of a TYPE II system. The problems associated 
with the pitch control pilot development stemmed from the presence of the high frequency 
complex poles (see EQ(IV.B.l.-2). These poles are quite close to the imaginary axis and 

due to the destabilizing distortions of the delay Root-Locus, can be easily shifted into 
unstable conditions. Closing the loop near w n = 3 rads/sec introduced large resonant 

reactions in the vehicle's lateral-directional body plane (X-Y body plane), due to torque 
reactive disturbances. To reduce this type of behavior, the closed loop dominant poles were 
placed near w n = 1.8 rads/sec and £ = 0.65. The pitch control pilot transfer function was 

given by: 


r pilot , , K pT (S + 0.8) e 0 2S 
U THETA ( ' S '' S(S + 10) 

The gain profile is listed in Table V. A. 1.-2. 

Additional Pilot Mechanisms for the Black Hawk 


(V. A. 1.-3) 


Tables V. A. 1.-3 through V. A. 1.-7 list the remaining single variable pilot mechanisms 
for the Black Hawk, that were developed in [18,19]. The models are based on the form: 
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.BHAWKf n K (S + A) e-°- 2S 
'pilots^ (S + B)(S + 10) 


(V. A. 1.-4) 


V.A.2. - Selection of Multi-Variable Configurations 

The maneuvering characteristics associated with flight control objectives will determine 
the configuration of the insertion set. The insertion set configuration is based on the multi- 
variable control structure shown in Figure III.B.3.-1 (for the Harrier, the Black Hawk 
provides only four cockpit control mechanisms [18]). Some typical flight control 
maneuvers that are considered in [18,20] are listed below. 

1. Pitch Reorientation 

2. Velocity Translation 

3. Altitude Translation 

4. Small Scale Heading Modification 

5. Altitude Rate Maneuver 

6. Flat Turn 

7. Coordinated Turn 

These maneuvers can be executed by the control configuration shown in Figure III.B.3.-2. 
The application of this method to this class of control maneuvers is straight forward. The 
insertion sets and their associated configurations for the above maneuvers are discussed in 
[18,20], and are summarized in Tables V.A.2.-1 and V.A.2.-2. 
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V.A.3. - Results of Static Pilot Insertion 


The results of pilot insertion can be best illustrated by considering several examples of 
multi-variable McRuer-Krendel pilot configurations actively participating in the control 
structure of the non-linear simulation program. Further results of static pilot insertion can 
be seen in [18,19,20]. 

Harrier : Pitch Reorientation 

The pitch reorientation maneuver is designed to provide pitch angle translation while 
maintaining level flight characteristics in the lateral, altitude, and velocity components. The 
multi- variable pilot configuration, for this maneuver, is given in Table V.A.2.-1. The pitch 
angle rotation will cause a redirection of the primary thrust vector, which will directly 
influence the altitude and velocity components. The velocity and altitude component pilot 
mechanisms will therefore be called upon to provide the necessary compensation. The 
command sequence issued to the pilot is of the form of a step in pitch angle while the 
remaining variables are held at their trimmed values. A 10 degree step maneuver in pitch 
from a 10 knot trimmed flight at 100 feet is driven by the control sequence shown in Table 
V. A. 3.-1. 

The command sequence was injected at 1 second into the simulation run. Plot V.A.3.-1 
shows the execution of the 10 degree set in vehicle pitch angle. This closed loop pitch angle 
response indicates a reasonably close match to the closed loop poles placed in the Root- 
Locus design. Plot V.A.3.-2 shows the pitch control pilots manipulation of the longitudinal 
stick. The change in pitch angle causes the primary thrust vector to be reoriented, creating a 
decelerating disturbance in the forward velocity. Plot V.A.3.-3 shows the forward velocity 
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response due to the pitch maneuver. Plot V.A.3.-4 shows the velocity pilot's compensative 
reactions to the velocity disturbance. Plot V.A.3.-3 shows that the pilot successfully 
reoriented the thrust vector and restored the commanded forward velocity. It is interesting 
to note that the final, steady state nozzle redirection is approximately 10 degrees (the 
commanded pitch angle). Plot V.A.3.-5 shows the altitude reaction due to the pitch 
maneuver. Plot V.A.3.-6 shows the pilot's manipulation of the throttle to handle the 
altitude disturbance. These plots indicate that the velocity pilot's reactions were fast enough 
to maintain proper thrust vector orientation with respect to the inertial frame (i.e. the 
velocity pilot was able to keep the thrust vector pointed towards the ground and therefore 
suffered little lift component degradation). The symmetry associated with this maneuver 
created insignificant disturbances in the lateral-direction modes. 

Harrier : Velocity Translation 

The velocity translations maneuver is designed to modify the forward velocity while 
maintaining level flight at a constant altitude. The multi-variable pilot configuration, for this 
maneuver, is given in Table V.A.2.-1. This maneuver uses the redirection of the primary 
thrust vector to modify the forward thrust. Reorienting the primary thrust vector tends to 
disturb the flight characteristics of the pitch and lift components. The command sequence 
issued to the pilot is of the form of a step in forward velocity while the remaining variables 
arc held at their trimmed values. A 5 knot step reduction in forward velocity from a 25 knot 
trimmed flight at 100 feet is commanded by the control sequence shown in Table V.A.3.-2. 

This command sequence was injected at 1 second into the simulation run. Plot V.A.3.- 
7 shows the execution of the 5 knot reduction in velocity. Plot V.A.3.-8 shows the velocity 
control pilot's reaction to the application of the command. Plot V.A.3.-7 's response 
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characteristic is very sluggish. This is due to the reduction in closed loop bandwidth via 
forward path gain reduction to reduce the pilot's control deflection. The initial pilot 
parameter selection created a response characteristics that tended to overdrive the nozzle 
angle (0j > 98.5°) during decelerating maneuvers. The velocity pilot parameters were 

adjusted to prevent the control mechanism overdrive, see Plot V.A.3.-8. Plot V.A.3.-9 
shows the pitch angle reaction to the thrust vector redirection during the execution of the 
velocity control maneuver. Plot V.A.3.-10 shows the pitch control pilot's compensative 
response due to the "nose down" effects of the velocity translation. The steady state error 
of Plot V.A.3.-9 is typical for the TYPE I system characteristic of the piloted pitch control 
loop. Plot V. A. 3.-1 1 shows a very small altitude response due to the velocity change. Plot 
V.A.3.-12 shows the small scale adjustments to the throttle by the pilot mechanism. 

Black Hawk : Coordinated Turn 

The coordinated turn maneuver provides a heading change through the execution of roll 
operations. The coordinated turn is executed be performing a roll maneuver while 
minimizing sideslip and maintaining level flight (pitch angle only) and constant altitude. 
The multi-variable pilot configuration, for this maneuver, is given in Table V.A.2.-2. It is 
interesting to note that the tail rotor collective pedal pilot mechanism has had it's visual 
feedback redirected from the heading to the sideslip indicator (slip ball). Heading control 
must therefore be supplied by the commanding process (i.e. the navigation/guidance 
process that is issuing the command). The roll angle modification will tend to cause a loss 
of the lift components due to the redirection of the main rotor thrust vector. The command 
sequence issued to the pilot is of the form of a step in roll angle, while the remaining 
variables are held at their trimmed values. A 20 degree bank turn from an 80 knot trimmed 
flight at 200 feet is driven by the control sequence shown in Table V.A.3.-3. 
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The command sequence was injected at the initiation of the simulation run. Plot V.A.3.- 
13 shows a smooth roll angle response to the command. The heading response, shown in 
Plot V.A.3.-14, indicates the execution of the turn. The transient heading response is due 
to the settling of the roll angle and the turn coordination operations. Plot V.A.3.-15 shows 
the suppression of the sideslip angle and indicates that turn coordination has been achieved. 
Plot V.A.3.-16 shows the altitude response during the bank turn and shows the 
compensatory effects of the altitude control pilot. Plot V.A.3.-17 shows the compensatory 
effects of the pilot mechanism on the pitch angle. 

V.A.4, - Conclusions on the McRuer-Krendel Models 

The McRuer-Krendel pilot models have shown their ability to operate the V/STOL 
aircraft over a wide range of flight control maneuvers. Control and regulation activities 
have been shown to adequately achieve the control objectives. The multi-variable 
configurations provide for a wide range of possible flight scenarios. The fixed structure of 
the multi-variable configurations does, however, create certain difficulties when the flight 
control objectives require the use of multiple configurations. The switch between 
configurations has been examined in [19]. Finally, the sluggish response characteristics 
due to the reduction of closed loop bandwidths to accommodate various vehicle constraints 
tends to limit successful prediction of human behavior. 
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V.B. - Optimal Control Model Pilots 


The concepts involved in the insertion of the OCM are based on the selection of two 
parameter sets: 1) Control task related, 2) Pilot related. The control task parameters 
describe the vehicle under control and the control objectives associated with the flight 
maneuvers that the pilot is required to perform. The pilot parameters describe the basic 
human response characteristics and inherent limitations. This section will illustrate the 
design and insertion of the OCM within the control loops of the Harrier and will present the 
results of various flight control tasks along with a comparison of the response 
characteristics of the OCM to some actual pilot flight data. 

V.B.l. - Description of the Control Tasks and Display Configurations 

The flight control tasks that have been considered for the testing of the OCM are 
classical precision hovering maneuvers that are performed outside of the ground effects 
region. Two primary maneuvers have been utilized to analyze the OCM: 1) Vertical tracking 
maneuver, 2) Lateral tracking maneuver. Each maneuver is complex in nature and exploits 
various aspects of the pilots control characteristics. The vertical hover maneuver consists of 
traversing between and positioning/aligning the vehicle at two vertically spaced targets, as 
shown in Figure V.B.l.-l. The targets, described in [25], are placed at 40 feet and 80 feet 
above the runway surface, and thus are outside of the vehicles ground effects region. The 
control objective is to hold at one target for a period of time, traverse to the other target at a 
constant rate, hold at that target for a period of time and then return to the initial target. 
Alignment with the targets is considered to be a positioning/disturbance regulation 
operation, while motion between the target is considered to be a rate control/tracking 
maneuver [25]. The Lateral tracking maneuver is similar to the Vertical task with the 
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exception that the targets are separated horizontally, as shown in Figure V.B.1.-2. The 
targets are placed 40 feet apart at an altitude of 50 feet. 

The vehicle that is used in this study is the Harrier II AV-8B. The cockpit controls 
available to the pilot, for the maneuvers in question, are: 1) Longitudinal stick (right hand), 
2) Lateral stick (right hand), 3) Throttle (left hand), 4) Rudder pedals (feet). The 
complexity of the maneuvers and their basic control structure, requires that the pilot 
maintain a fixed nozzle angle, thus forward velocity and longitudinal positioning control are 
indirectly handled through pitch angle control. 

The displays that the pilot monitors are assumed to be the targets (external visual cues). 
The targets have been designed [25] to supply the pilot with the following information: 
Pitch, Roll, and Yaw angles, Longitudinal, Lateral, and Vertical positions. The pilot can 
therefore derive the output vector shown in Table IV.A.3.-1. The state space model and the 
above output vector, derived in Appendix B, are based on this display assumption. 

V.B.2. - Overview of the Piloted Flight Data and Analysis 

The precision hover maneuvers, described above, where used in a simulation fidelity 
study by [25]. The time domain flight test data of the actual pilot activities and the vehicle 
responses from this study were supplied by NASA-Lewis. This information provided a 
valuable tool in the design and insertion of the OCM pilot. The flight data was reviewed 
and the basic flight maneuver trajectories were extracted. This trajectory information is the 
basis of the command sequences that are injected into the OCM. 
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V.B.3. - Description of Pilot and Task Parameter Selection 

The OCM task parameters were selected from an analysis of the fundamental structure 
of the control tasks in the piloted flight data. The limited information available from the 
flight data required that most task parameters were selected subjectively. Thus, a simplified 
performance index was chosen to represent the task of minimizing the hover attitude and 
position errors indicated by the external target. 



0 


(V.B.3.-2) 

The maximum desirable values of the position errors (XX M , YY M , ZZ M ), were chosen to 

correspond to the relative target size and reflects the specifications associated with precision 
hovering. The values of the attitude angle cost parameters were chosen subjectively by 
considering a 10 degree deviation as being large. The maximum values of the cockpit 
control mechanisms were selected to be approximately 0.25 of the total control travel, as 
suggested in [26]. Table V.B.3.-1 and V.B.3.-2 list the values of the cost function weights 
associated with the vertical and lateral flight maneuvers, respectively. 
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The vertical maneuver parameters were originally designed for all maneuvers. During 
tests (discussed later sections), the lateral operations showed unstable response 
characteristics in the pitch and yaw (heading) components when using the values of Table 
V.B.3.-1. Tightening the allowable deviations in pitch and yaw (compare Tables V.B.3.-2 
to V.B.3.-1), improved these response modes. The response characteristics in the yaw 
angle are expected because of the inherent couplings in the lateral-directional components. 
The responses of the pitch angle were more pronounced than were expected. These 
problems appear to stem from the limitations of the OCM’s internal reference model, since 
the optimal control gains are based on the internal reference model (it's knowledge of the 
system under control). Another possible reason for the differences of Tables V.B.3.-1 and 
V.B.3.-2 can be reflection of the OCM's sensitivity to the inherent differences between the 
two maneuvers. This may be due to a reconfiguration of the manner in which the pilot must 
obtain information during the maneuvers. 

The pilot parameters were chosen according to some typical values. The weightings of 
the control rates were selected to provide a neuro-muscular lag time constant, T N , near 0. 15 

seconds for each control mechanism. The magnitude of the motor noise sources for each 
control is shown in Table V.B.3.-3. These correspond to an approximate -15dB signal to 
noise ratio (S/N). The magnitude of the observation noise was chosen to correspond to an 
approximate -lOdB (S/N). These values create a remnant that is considered reasonable for a 
multi-axis hovering control task [10]. The pure time delay was selected to be 0.2 seconds. 

V.C.4. - Results and Comparisons of OCM Pilot Insertion 

The OCM pilot was inserted into the nonlinear simulation environment of the Harrier. 
The structure of the OCM required that the vehicle trajectory associated with the control 
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task be specified in position/attitude and velocity/rate. This type of trajectory is generated 
by the velocity/rate driven system shown in Figure V.B.4.-1, for the pitch angle 
components. This generator is replicated for the other state variables. The low pass filter 
provides bandwidth limiting of the forcing function as described in [2]. The time constant 
of this filter was selected to be 2.5 seconds. 

An analysis of the vertical and lateral flight data showed a variety of pre-test procedures 
used by the pilot to approach and align the vehicle with the target. In both cases, the vehicle 
was at rest on the ground away from the targets. The approach and alignment operations 
consisted of a vertical takeoff, followed by a ground translation and a flat rotation (yawing 
maneuver). During initial attempts at constructing the command sequence for these 
maneuvers for the OCM, trim discrepancies in the non-linear simulation program created 
problems. These primarily stemmed from problems associated with simulating the vehicle 
on the ground. To overcome these problems, the OCM tests were initialized near the targets 
with the vehicle already airborne. This allowed closer matches of trimmed values and 
simplified the test command sequences. 

Vertical Tracking Maneuver 

The vertical tracking maneuver consisted of an initial yawing rotation (to simulate the 
final stages of the pilot's alignment with the targets), followed by a holding period and then 
by the cyclic execution of the vertical operations. The cost function weightings of Table 
V.B.3.-1 were used to select the control aspects of the OCM. The vehicle was trimmed to 
the values shown in Table V.B.4.-1. The vertical rate and yaw rate command sequences of 
Figure V.B.4.-2 were applied to the OCM pilot during a 125 second simulation run. This 
command sequence requires the pilot to rotate the vehicle (flat tum - approximately 60°) to 
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acquire the upper target. The pilot holds at the target (for approximately 10 seconds) then 
proceeds to the lower target at a constant rate (2.5 feet/sec). The pilot is then repeatedly 
commanded to proceed back and forth between targets. 

Plot V.B.4.-1 shows a heading angle comparison of the OCM (solid line) and the 
piloted flight data (dashed line). This type of comparison strategy will be used throughout 
the following discussion. Transients associated with the target alignment phases are not 
present for times greater than 35 seconds (TIME>35 seconds). Plot V.B.4.-2 shows a 
comparison of the rudder pedal activity. The "fuzziness" of the OCM response is due to the 
noise model of the controller remnant. Rudder pedal operations differ significantly. This 
has been attributed to the limitations in the lateral dynamics of the OCM's internal reference 
model. The flight data tends to suggest that the heading angle regulation is not particularly 
critical to this maneuver. It appears that this type of activity may stem from a threshold 
regulation process where no corrective actions are taken until the error crests a certain level. 
The present version of the OCM does not account for this type of behavior. 

Plots V.B.4.-3 and V.B.4.-4 compare the altitude and vertical rate responses, 
respectively. These plots show similar trends in command trajectory following. Target 
alignment operations can be seen during the initial phases of the run (TIME<35 sec). The 
vertical rate responses of Plot V.B.4.-4 show that the OCM has similar behavior during the 
execution of the vertical maneuvers (TIME>35 sec), but does not capture the higher 
frequency content of the pilot data. This has been attributed to the low order altitude 
component representation of the of OCM's internal reference model, and possibly the lack 
of auditory feedback. 
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Plot V.B.4.-5 shows a very dissimilar comparison of throttle operations. Plot V.B.4.-6 
shows a comparison of the engine speeds. During the execution of the vertical maneuvers 
(TIME>35 sec), the engine speeds show strikingly similar response characteristics from 
very dissimilar throttle operations. This discrepancy may be due to infidelities in the 
simulation environment, particularly in the throttle linkages and servo-systems. 

Plot V.B.4.-7 compares the pitch angle responses. This plot shows very similar vehicle 
orientations during the vertical maneuvering. This is an important response characteristic 
due to the fixed nozzle angle which couples the pitch angle to the forward velocity and 
position. Plot V.B.4.-8 shows the comparison of the longitudinal stick activity. The offset 
in stick operations is due to the trimmed value of the piloted data being —1.35 inches. As in 
the case of the rudder pedal, the OCM does not predict the higher frequency operations of 
the pilot. 

Plot V.B.4.-9 shows a comparison of the airspeed responses. The OCM pilot shows 
very similar characteristics in both magnitude and frequency content during the vertical 
maneuvers (TIME>35 sec). Plot V.B.4.-10 compares the longitudinal position. The offset 
is due to the initial target approach operations by the pilot. This plot shows similar but out- 
of-phase positioning response. 

Plots V.B.4.-11 and V.B.4.-12 show comparisons of the roll angle responses and 
lateral stick deflections, respectively. In both plots, the OCM fails to match the frequency 
content of the actual pilot. 
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Lateral Tracking Maneuver 


The lateral tracking maneuver was simplified by not including the initial yawing 
operations of the target alignment. The vehicle was trimmed to the values shown in Table 
V.B.4.-2. The lateral velocity command sequence of Figure V.B.4.-3 was applied to the 
OCM pilot during a 100 second simulation run. 

Plot V.B.4.-13 shows a comparison of the lateral position responses. This plot 
indicates a good trajectory following by the OCM. Plot V.B.4.-14 compares the lateral 
velocity responses. As in the vertical maneuvers, the OCM’s response characteristics do 
not capture the higher frequency content of the pilot data. This is illustrated in the 
regulatory operations shown in Plots V.B.4.-15 - V.B.4.-20. 

V.C.5. - Conclusions on the OCM pilots 

The OCM pilots have shown their abilities of providing adequate multi-variable flight 
control of the V/STOL research aircraft. The OCM is not structurally limited in it's 
maneuvering characteristics and can therefore be applied to a wide range of flight control 
objectives. Differences between the cost function weightings of the vertical and lateral 
tracking maneuvers illustrate the pilot's alteration of his approach to the individual control 
problems (possible modification of feedback to handle the characteristics of a specific 
maneuver). The OCM showed an excellent match in the longitudinal components of the 
vertical maneuver. This appears to stem from a reasonable model of the aircraft's 
longitudinal dynamics (since lateral excitation is minimal) and subjective assumptions in the 
cost function weight selection. 
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The reasons behind the OCM's inability to predict the higher frequency components of 
the pilot data are not clear. A possible reason stems from the manner in which the pilot 
obtains information. The limited information available on the targets and actual flight 
control objectives of the piloted flight data required the use of an intuitive selection of 
displays. A close examination of the pilot’s control activities showed a higher degree of 
sensitivity to the angular rates than to the angular position. Attempts at configuring the 
OCM to be more sensitive to angular rates did not resolve this issue. The addition of a 
display scanning algorithm induced a cyclic feature in the OCM's responses that was 
similar to that of the pilot data. The presence of scanning behavior in the pilot flight data 
can not be substantiated and thus the OCM scanning model is only marginally permissible. 
Another possible reason is the OCM's remnant models. In some regulation operations [10], 
the pilot tends to "battle" his own remnant more than the external disturbance. Modifying 
the OCM's remnant models did not significandy alter the frequency content. Finally, it is 
possible that the OCM's reliance on the internal reference model of the vehicle dynamics 
creates this phenomena. Limitations in the model may restrict the manner in which the 
OCM applies it's control efforts. The lack of deadband considerations in the model (i.e. the 
0.1 inch deadband in the servo-channels of the longitudinal and lateral stick linkages) leads 
to the suppression of OCM control mechanism deflections that are based on small errors. 
Possible over-accuracy in the model may provide the OCM with such a good understanding 
of the vehicle under control that all compensative actions are near optimal. This would 
account for the OCM's limited control deflection activities in the regulation modes. 

The OCM is very sensitive to changes in cost function weightings and muscular system 
time constants. In some cases (pitch and roll components), maximum deviation values of 


61 


20-30% larger produced unstable results. This may be due to the nonlinearities in the 
simulation environment. This tends to complicate the subjective selection of cost functions. 
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VI. - SUMMARY AND CONCLUSIONS 


This report has illustrated the application of "Paper Pilots" within the flight control 
loops of V/STOL research aircraft. Two types of human dynamics models have been 
considered: 1) McRuer-Krendel - single variable transfer function, 2) Optimal Control 
Model - multi-variable approach based on optimal control and stochastic estimation theory. 
Descriptions of the models and discussions of their inherent limitations in predicting human 
behavior have been provided. Design strategies and methods of inserting the pilots within 
the control loops have been discussed and illustrated for two V/STOL research aircraft: 1) 
Sikorski Black Hawk UH-60A - a high performance helicopter, 2) McDonnell Douglas 
Harrier II AV-8B - a thrust vectored jet fighter. Results of simulated pilot insertion have 
been analyzed and compared to the control activities of actual human pilots performing 
similar control objectives. 

The "Paper Pilots" have shown their abilities to successfully "fly" the V/STOL aircraft 
that have been considered in this research. The simulated pilots provide a stabilized control 
over the vehicle and respond to control objectives in a manner similar to human pilots. The 
response characteristics of the pilot models are, however, very general in nature. This 
appears to stem from the "Paper Pilots" inability to completely simulate the human s 
complex manner of obtaining and processing information, and applying the controls to the 
vehicle. In addition, it may be possible that the pilots have an over-detailed description of 
the vehicle under control and thus exert a nearly ideal control response. 
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Of the two human dynamics models that have been considered, the OCM appears to 
provide the better simulation of human activities. It's multi-variable structure and state 
variable framework provide a simple yet effective method for describing the pilot's abilities 
and the control objectives that the pilot is confronted with. The maneuvers that have been 
considered for the OCM are in the low speed/powered-lift region of the flight envelope. To 
maneuver in other regions of the flight envelope, the control objectives (i.e. the vehicle 
model and control strategies) must be altered because of the changes in vehicle dynamics. 
The use of envelope based gain tables provides a possible solution but the tables must be 
constructed in an off-line approach (not a simple task) and do not allow for alternative flight 
configurations. In addition, transition/interpolation between flight regions must be well 
defined to avoid control transients. 

Improvements are needed to provide the OCM with a more complete description of pilot 
activities and to simplify the OCM pilot design process. Adaptive procedures have been 
successfully applied to the single variable McRuer-Krendel model. This strategy appears to 
have a good prospects for use in the OCM based pilots. The OCM is, however, a much 
more detailed and complex system. Adaptive procedures for the OCM will require a 
different parameter estimation scheme (because of it's high order) and will also need a 
different optimization scheme if muscular system time constants are to remain the same. 
Initial research on this subject (not discussed in this report) indicated favorable results. Due 
to time constraints, this work was not fully established. A further investigation in these 
areas is suggested as the focus of future research efforts. 


In conclusion, the "paper pilot" provides a safe, convenient, and seemingly effective 
method for introducing human pilot response characteristics during the design process 
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without resorting to the use of manned, fixed-base simulators. Construction and application 
of the pilot models are relatively straight forward and the pilots can be configured to 
achieve a large number of flight configurations and control objectives. 


65 



APPENDIX A 


DESCRIPTION OF THE OCM 

This appendix presents a mathematical overview of the Optimal Control Model of 
Human Dynamics that is developed in [4]. A continuous time description of the OCM is 
discussed.and a discrete time representation is then developed for use in the computer 
simulation environment. 


A.l. Mathematical Overview of the OCM 

This section provides a general overview of the mathematical concepts involved with 
the OCM. Figure A.l.-l shows a block diagram of the internal structure of the OCM. For 
the purpose of definition, the dynamics of the system under control will include the 
dynamics of the actuation, sensory subsystem, SAS, and any other on-board control 
systems. The overall system under control will be represented by a set of linearized 
equations of motion. 

x(t) = Ax(t) + Bu(t) + w(t) (A.l.-l) 

y(t) = C x(t) (A. 1.-2) 

The system dimension is assumed to be "n", the number of inputs (i.e. cockpit control 
mechanisms) is assumed to be "m" and the number of outputs (i.e. visually and extracted 
displays) is assumed to be "1". The variables and parameters are defined as follows: 
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x(t) = "n" dimensional system state vector 


u(t) = vector of "m" control inputs (in this case, the cockpit control mechanisms) 

y(t) = vector of "1" system outputs (linear combination of the system states as 
perceived and deduced by the human from the displayed information) 

w(t) = vector of "n" external disturbances that are independent, zero mean, 
white, gaussian noise sources, where: 

E{w(t),w T (a)) = Q w 5(t-a) (A. 1.-3) 

A = (n,n) system matrix 
B = (n,m) input distribution matrix 
C = (l,n) output measurement matrix 


It is assumed that the human pilot/operator maintains an internal model of the system 
under control to base his control and estimates on. 


x(t) = Fx(t) + G u(t) (A. 1.-4) 

y(t) = H x(t) (A. 1.-5) 

The parameters are defined as follows: 

F = Human's perception of the system matrix 
G = Human's perception of the input distribution matrix 
H = C 
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It is assumed that the control task is adequately reflected in the human's choice of the 
best control input, "u*(t)". In addition, the human's choice is also based on his inherent 
knowledge of his neuro-muscular limitations. The optimal control input "u*(t)" minimizes, 
in steady state, the cost function given by: 

T 

J(u) = T I ™ E | |J[ x(t) T Qx(t) + u(t) T Ru(t) + u(t) T Su(t)]dt ^ J (A.I.- 

6 ) 


where 

a< t 

Q and R are positive semi-definite 
S is positive definite 

The formulation of EQ(A.l.-6) does not directly include the neuro-muscular dynamics 
of EQ(IL-l), but instead provides a cost on the control rate in J(u). The inclusion of a 
control rate term results in a first-order lag being introduced in the optimal controller. This 
term is utilized to indirectly model the physiological limitations of the rate at which a human 
can perform a control action due to the neuro-muscular/motor dynamics. 

The solution to this optimization problem is obtained by defining an augmented system. 
i 0 (t) = F O *o(0 + G 0 u(t) + w 0 (A. 1.-7) 


where 
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(A. 1 .-8a) 



x(t) 

H(t) 


T 

G' 

.0 

0 . 



w 0 


w(t) 

0 


An optimal control rate is generated by the linear feedback law: 


u*(t) = - Kq X 0 (t) = - [ K° 0 KjJ 


x(t) 

u(t) 


From EQ(II.A.2.A.-3) and Figure H.A.2.A.-2 we have: 


T n i*(t) + U*(t) = u c (t) + V uu (t) 
u c (t) = - K* x(t) 

thus 

u*(t) = - u*(t) + T n i U c (t) + T n ’ V uu (t) 

or 

u*(t) = - T n * U*(t) - T n i K* x(t) + T-J V uu (t) 
which results in 


(A. l.-8b) 

(A. 1.-9) 

(A.l.-lOa) 

(A.l.-lOb) 

(A.l.-lOc) 

(A.l.-lOd) 
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(A.l.-lOe) 


u*(0 = * Kq X 0 (t) = -[Kj KjJ 


x(t) 

u*(t) 


+ Vuu(t) 


where 


K o 




= T-J K* 


Tn 1 


(A.l.-lOf) 

(A.l.-lOg) 


The feedback gain matrix, Kq, is obtained from: 


Ko = G 0 T S-> P 0 


(A.1.-11) 


Pq is the unique positive definite solution of the "n+m" dimensional matrix Ricatti 
Equation. 


F 0 ' Po + Po F o + Qo - PoGoS 'Go T Po = 0 


(A. 1.-12) 


where 


Qo = 


Q 0 
.0 RJ 


(A. 1.-13) 


The key to implementing this control strategy is in the selection of the control rate cost 
weightings, S, to obtain: 


K« - [ K °0 T n] 


(A. 1.-14) 
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The Kalman Filter estimates the delayed system state from the observation of the 
delayed, noisy system outputs and an inherent knowledge of the delayed "command” 
control, u c (t-d). The Kalman Filter uses an alternate augmented system to derive it's 

estimates. 


Xj(t) = Fj Xj(t) + GjUgCt-d) + G n W, 


(A. 1.-15) 


where 



x(t-d) 


w(t-d) 

Xj(t-d) = 


W! = 



u(t-d) 


V^t-d) 


(A.l.-16a) 



T 

G ' 


'O' 


'I 

O' 

Fi = 

.0 

- T 'n_ 

G, = 

Jn_ 

Gn = 

.0 

T n_ 


(A.l.-16b) 


Hj ■ [ H 0 ] 


(A.1.-16C) 


It is important to note that the estimated state of EQ(A. 1.-15) utilizes the delayed desired 
control input, u c (t-d), as the driven deterministic input. The Kalman Filter generates the 

delayed estimate of the delayed state via: 

Xj(t-d) = Fj Xj(t-d) + GjUcCt-d) + K t [ y p (t) - H t Xj(t-d) J (A. 1.-17) 
The output error is weighted by the Kalman gains, K t , that arc generated by: 

K, = Pj Hj 7 Ry 1 (A. 1.-18) 
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where the error covarience matrix, Pj, satisfies 

F, T P, + P, F, + Q, - P, H, Ry H, t P, = 0 


Q. 



O' 


where 

= E { T„> V„(t) , V uu (t) T T‘ n 1 T ] = T'jj Q u T„' t 


(A. 1.-19) 


(A.l.-20a) 


(A.l.-20b) 


The least-squared error predictor generates a time advanced estimate of the present time 
state from the delayed estimated state of EQ(A.1.-17). The prediction is based on the 
internal reference model: 


C(t) = Fj C(t) + G lUc (t) 


where 


C(t) = 


x(t) 

u c (t) 


The time update of the predicted state is generated by: 

Xj(t) = f(t) + K 2 [ Xj(t-d) - f(t-d)] 

The delayed state error weighting, K 2 , is given by: 


(A. 1.-21) 


(A. 1.-22) 


(A. 1.-23) 
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(A. 1.-24) 


which is the state transition matrix for the time advancement. 

The fundamental modules of the OCM, that have been defined above, are summarized 
by Figure II.B.-2. 


A.2. Discrete Time Representation of the OCM 

The continuous time OCM described above was transformed to a discrete time 
equivalence for insertion into the simulation environment. The following sections describe 
the transformation technique and the resulting discrete time model. It is important to note 
that the Kalman Filter and the Optimal Control gain calculations will be based on time 
varying solutions. 

Discrete Time System Model 

The discrete time representation of the system’s involved in the OCM are based on the 
Zero-Order- Hold equivalence transform (ZOH) described in [27]. The difference equation 
representation of the perceived system’s state equation, EQ(A.l.-4), is given by 

x k+1 = <&x k + Tu k + T w w k (A. 2.-1) 

and the measurement equation, EQ(A.l.-5), is given by 
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y k = H x k + Hy V k 


(A.2.-2) 


The discrete time "n" system states, "m" system inputs, and "1" system outputs are 


represented by 

x k = [x k ,x k ,...,x k ] T (A. 2. -3) 

u k = [u^,u^,...,u^] T (A.2.-4) 

yk = [yk-yk ykl T (A. 2 .- 5 ) 

The discrete time process noise model and covarience matrix are given by 

w k = [w k ,w k ,...,w k ] T (A. 2. -6) 

QZ = E [ w k , w k T ] = r w Q w r w T (A.2.-7) 

The discrete time observation noise model and covarience matrix are given by 

v k = [v},v2,...,v y n ] T (A. 2.-8) 

= E [ V£ , Vf ] - Ry (A.2.-9) 

The ZOH discrete time equivalence transformations are given by 


(A.2.-10a) 


74 



(A. 2. -10b) 


T s 

r = Je FTs dCT G 
T 

A S 

r w = Je FT *da (A. 2. -10c) 


where 

T s = Sampling period 

Pure Delay Model 

Information processing and neuro-motor signal delays are represented by the pure 
delays shown in Figure A.l.-l. The delays are implemented by a sliding window FIFO 
buffer of length D. 

D = MOD(^) (A. 2.-1 1) 


where 

D = Number of sampling periods 
d = Continuous time delay 
T s = Sampling period 

The FIFO buffer is driven by 

DO I - D. 1,-1 

X KJ = X K . W (A.2.-12) 

ENDDO 
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Neuro-muscular and Optimal Control Generator 


A recursive time varying solution of the discretized optimal control/gain generator 
described in [27] is utilized. The optimal control rate is defined by EQ(A.l.-lOe) and it's 
discrete time representation is given by 

u*k+i ” u k + r K 0 x k + T k i V“ (A. 2. 13) 

where the motor noise and it's covarience matrix are given by 

Q! = E[vs,vf] = r K .Q„r Kil T 

The optimal state feedback gains can be rewritten into the form of the recursive solution 


(A.2.-14a) 

(A.2.-14b) 


u k+l 



+ r KX 


(A.2.-15a) 


where 


K - [ K ” k <M] 


(A. 2. -15b) 


The feedback gain matrix is obtained from the recursive relations in terms of the augmented 
system 
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(A.2.-16) 


On 


o r 

0 0 . 


r o 


'O' 

. 1 . 


The recursive solution of the optimal control gains is given by: 


k» = (s k + r 0 T p»r 0 )'r 0 T p»«. 0 

(A. 2.- 17a) 

Pfc = Pk - P£ r o( S k + r 0 T Pk r o) lr o T ^k 

(A. 2. -17b) 

P» = <t 0 T P«<P 0 + Q» 

(A.2.-17C) 

The cost weightings on the states, controls, and control rates (Q,R,S) respectively, are 
introduced in the manner used in the continuous problem 

TQ °1 

QP = s k = S 

^ Lo Rj 

(A.2.-18) 


The control rate weightings, S k , are chosen such that, in steady state, the set of equations 
EQ(A. 1.-17) result in 


r- T s 

e' T n 


-1 


<l*„] = e - K i T s = e' T = T „ . 

*0 


T s 

e" T n 


e t Jj 1 _l 


(A.2.-19) 


The state feedback matrix directly results from the solutions of EQ(A. 1.-17). The diagonal 
elements represent the time constants of the individual neuro-muscular systems associated 
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with the separated cockpit control mechanisms. The motor noise gain matrix is obtained 
from 


r„i = I-<V (A.2.-20) 

The finalized state feedback matrix results from 

r d = r K^ Kok (A. 2. -21a) 


Discrete Time Kalman Filter 

A recursive time varying solution of the discrete time Kalman Filter described in [28] is 
utilized. The Kalman Filter is used to estimate the delayed system state from the 
observation of the delayed noisy system outputs 

yk = yk-D + v ki> (A. 2. -22) 

from a knowledge of the desired delayed control command 

A A 

= [ u l> u 5>”-> u m]k-D T = ‘ K -0 x k = ' *k (A. 2. -23) 

and an understanding of the motor and process noises and their respective covariences. The 
Kalman Filter bases it's estimates on the augmented system given by 
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F“A 

X 


X k-D ~ 


k-D 


\-D 


w 1 

W k-D 


r~ a -i 

^k-D 

A 

V u 

V k-D 


o x = 


(D r 

o % 


r, = 


- o 


(A. 2. -24a) 


(A. 2. -24b) 


The measurement update equations are 


X k 'n - X‘ + Kt[yJ - HX>] 


P‘ = PJ - K>HP> 


(A. 2. -25a) 
(A.2.-25b) 


and the time updates are given by 


x k !o = <*>> x k ' D + r, u t ; 


k-D 


P> - O, PiV + Qi 


(A. 2. -26a) 
(A.2.-26b) 


where 


<4 = 


Q 1 

o 


O' 

A 


Qk “ ^1 Qk ^ 


(A. 2. -26c) 


The time varying Kalman gains are generated by 


Kj. = P£H T (HP£H T + ) 1 


(A.2.-27) 


with the initial conditions 
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(A.2.-28a) 


Pi = 0 

A _ 

X k l D = ^ (A. 2. -29b) 

Discrete Time Predictor 


The discrete time, time advance predictor uses the internal reference model given by 


Cm - «i Ck + r i 


rA 



X k-D 

A 


U k-D 


The prediction matrix is given by 

A A 

xl - C k + K l [ X t ] D - 4-d 1 

where the state projection matrix is given by 



(A. 2. -30a) 


(A. 2. -30b) 


(A.2.-31) 


(A.2.-32) 
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c.z- 


APPENDIX B 


A LINEAR STATE SPACE MODEL 

This appendix presents a linear state space model of the Harrier AV-8B dynamics from 
trim, while operating in the low speed/powered lift region of the flight envelope. This 
model is directly utilized as the internal reference model in the OCM implementations. This 
appendix is separated into two sections: 1) Assumptions, simplifications and general 
derivation of the model, 2) Identification of model parameters from an analysis of the 
response characteristics of the non-linear simulation program. 

B.l. - Harrier AV-8B Model Development 

A generalized linear representation of the core longitudinal dynamics is given by: 
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(B.1.-6) 


Due to the low speed assumption, a variety of simplifications can be made and verified by 
an analysis of the responses of the simulation model program. The first series of 

assumptions are directed at the forward velocity components. A very general assumption 
that can be made is: X q ~ 0. This assumption is based on the relative unimportance of the 

angular rate on the forward velocity through the entire flight envelope [21,22]. The next 
simplification is: X w ~ 0. This assumption is based on aerodynamic symmetry and an 

analysis of the vehicle's forward velocity response characteristic due to an positive impulse 
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of the throttle while in a near hover as shown in Plot B.l.-l. An additional simplification 
of X T ~ 0, can be derived from Plot B.l.-l and an intuitive analysis of the low 

speed/powered-lift thrust vector's effect on the forward velocity as shown in Figure 
IV. A. 1.-2. The final forward velocity component assumption is: X e ~ 0. This 

simplification is based on the longitudinal stick's dominant effect on vehicle pitching 
motions, shown in Figure IV.A.1.-3. Plot B.1.-2 backs-up this assumption by illustrating 
that the longitudinal stick's effect on forward velocity is coupled to the redirection of the 

thrust vector associated with the pitching motion and is therefore not direct. Thus as shown 
in the thrust diagram of Figure IV.A.l.-l, the nozzle angle control, Xj, will dominate the 

control of the vehicle's forward velocity. It is important to note that this model does take 
into account the pitch angle's effect on the forward velocity, where typically, X 0 = g, the 

gravitational acceleration. 

The effects of the lateral system components on the response characteristics of the 
forward velocity where considered insignificant from a analysis of the responses of the 
simulation model program. Plot B.1.-3 shows the forward velocity response due to an 
impulse on the lateral stick. Plot B.1.-4 shows the response of the same state variable due 
to an impulse on the rudder pedals. The forward velocity reactions to the lateral stick and 
rudder pedal operations are similar in magnitude to those of the longitudinal stick and 
throttle. For this reason, the lateral components of the forward velocity are ignored. 

The next series of assumptions are directed at the rotational modes of the longitudinal 
dynamics. The first assumption is based on aerodynamic symmetry and the negligible 
pitching reactions due to vertical motion, thus ~0. Plot B.1.-5 illustrates the relatively 

small pitch rate response due to an impulse of the throttle. Although M w does have a small 
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effect on the pitching rate, it will not effect the hovering characteristics. An additional 
assumption that can be derived from Plot B.1.-5 is M T ~ 0. This term primarily describes 

the moment arm of the thrust vector on the vehicle's center of gravity. A final rotational 
component assumption is Mj ~ 0. This term describes the rotational effects of the nozzle 

angle due to the thrust vector's moment arm on the vehicle center of gravity. Plot B.1.-6 
shows pitch rate reaction due to an impulse of the nozzle angle. Again, this response 
characteristic will not have a significant effect on the low speed flight dynamics. The 
control mechanism assumptions of M T ~ Mj ~ 0, indicate that the longitudinal stick will 

dominate the control of the rotational dynamics due to the relatively large moment arm 

associated with the physical locations of the forward and aft RCS jet vents, as shown in 
Figure IV. A. 1.-3. The term M u is retained to provide a coupling between the forward 

velocity and longitudinal rotations. Although, low speeds are assumed, it will be shown 
later that this term and the forward velocity terms, X e and X u , provide the couplings that 

tend to generate the long term Phugoid responses. 

The pitch rate responses due to impulses of the lateral stick and rudder pedals are 
shown in Plot B.1.-7 and Plot B.1.-8, respectively. These reactions are similar in 
magnitude to those of the throttle component, Mj, and are therefore neglected. 


The final series of assumptions in the longitudinal dynamics are directed at the vertical 

velocity components. The first assumption is based on the relatively small contribution of 
the pitch rate on the vertical components [21,22], thus, Zq ~ 0. Plot B.1.-9 illustrates the 

vertical rate response due to an impulse on the longitudinal stick. This response 
characteristics reveals that the longitudinal stick's component, Z e ~ 0, is also a realistic 

assumption. The pitch angle component is small due to the small angle and low speed 
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assumptions, thus, Z 0 ~ 0. Another simplification due to the low speed assumption is Z u ~ 

0. This assumption would not hold true if the Harrier was a tilt-wing V/STOL aircraft. For 
those types of vehicles, Z, u approaches Z w . The final simplification is directed at the relation 

between the nozzle angle and the vertical components. Plot B.1.-10 shows the response of 
the vertical rate due to an impulse of the nozzle angle. An examination of Plot B.1.-10 and 

Figure IV.A.l.-l shows that perturbations about the large nozzle angle has only a small 
effect on the lift components, thus Z- ~ 0. The control mechanism assumptions of Zj - Z e ~ 

0, indicate that the throttle control will tend to dominate the vertical dynamics due to the 
large nozzle angle associated with low speed/powered-lift flight. 


The vertical rate reactions due to operations of the lateral stick and rudder pedals can be 
seen in Plot B.l.-l 1 and Plot B.1.-12, respectively. These responses are on the order of 
those associated with the longitudinal stick and nozzle angle. This suggests that the lateral 
components can be neglected. 


Incorporating the above assumptions within the longitudinal system model of EQ(B.l.- 
6) results in the following low speed/powered-lift longitudinal model. 
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(B.1.-7) 




L&rJ 


To better understand the longitudinal responses of the low speed/powered-lift region of the 
flight envelope and the correspondence of the vehicle model that has been developed, an 
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analysis of the vehicle's longitudinal dynamics will now be conducted. This analysis is 
based on the Laplacian approach used in [22]. An additional objective of this analysis is to 
obtain an understanding of the transfer function relationships to the above system model. 
This will be helpful in implementing parameter identification techniques. Before initiating 
the analysis it should be observed that the vertical components are strongly decoupled from 

the translational and rotational components. This is primarily due to the dominance of the 
powered-lift assumptions associated with the vertical components and being negligibly 

small while in low speed flight. 

The rotational and translational components of the longitudinal dynamics can be 
approximated by the following system. 
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(B. 1.-8) 


The characteristic equation of this reduced system is given by: 


\ong = detCSI - A) = det 


s - x u -X e 0 
Os - 1 

-Mu 0 S - H _ 


which results in 


\ong = S(S - M q )(S - X u ) - X 0 Mu 


(B.1.-9) 


(B.l.-lOa) 
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or 

a long = ' (X„ + M q )S 2 + MqXyS - (B. 1 1 Ob) 

EQ(B.l.-lOb) is the classical longitudinal hovering cubic [22], From Plot B. 1.-13, one can 
observe both the short period and long period dynamics of the pitch rate's response. This 
type of response supports a characteristic equation of the form: 

- (S + T sp )(S 2 + 2w ]p S, p S + w, p 2 ) (B. 1 . ■ 11) 

where the parameters are given by: 

= The pole associated with the time constant of the short period response 
w lp = The natural frequency of the long period response 

5] p = The damping ratio of the long period response 
Expanding EQ(B.1.-1 1) to the form of EQ(B.l.-lOb) results in: 

S 3 + S 2 (2w lp S, p + T„) 4- S(w lp 2 + 2w lp 5, p V + T sp w lp 2 
Relating the terms of EQ(B.l.-lOb) and EQ(B. 1.-12) we have: 

-(X u + M q ) = 2w lp 5 lp + T sp 

M qXu = W lp 2 + 2w lp 5 lp T sp 

-X e M u = T sp w lp 2 
From [16,17,19], we can assume: 


(B. 1.-12) 

(B.l.-13a) 

(B.l.-13b) 

(B.l.-13c) 
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(B.l.-14a) 


T_ » w, 

sp lp 


w lp 2 is small (B.l.-14b) 

Applying these assumptions to EQ(B. 1.-13) results in: 

M^-Tsp (B.l.-15a) 

X u ~ "2w lp 5 lp (B . 1 . - 1 5b) 

Thus we can see that the long period modes are due to the coupling between the 

translational and rotational components. The short period modes are purely due to the 
rotation of the vehicle about it's center of gravity. X 0 transmits the short and long period 

rotational perturbations to the translational motion. My primarily couples the long term 

translational perturbations to the rotational dynamics. It can be shown, through a more 
detailed analysis [23], that M u < 0 causes divergence in the phugoid modes. This may be 

the case for the Harrier configuration that is used in this study. Plot B. 1.-14 shows a 
diverging phugoid mode in the pitch rate. This type of unstable phugoid oscillation appears 
to be typical for hovering vehicles when out of ground effect [22]. 


The resulting low speed/powered-lift longitudinal model for pilot frequencies is given 
by: 
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(B. 1.-16) 



87 



The characteristic equation of this system is given by: 


A long - S 2 (S - M q )(S - Z w ) 


(B.l.-17a) 


or 


^ LONG 


S 4 - S 3 (M q +Z w ) + S 2 M q Z w 


(B.l.-17b) 


A generalized representation of the lateral-directional dynamics is given by: 
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In a manner similar to the longitudinal model, a variety of simplifications can be made and 
verified by an analysis of the simulation model responses. The first series of assumptions 
are directed at the lateral velocity components. The coupling of the roll rate to the lateral 
velocity, Y p , can be neglected due to it's general unimportance throughout the flight 

envelope [21,22]. The primary coupling of the roll rotational dynamics is provided by the 
roll angle component, Y^. This is because of the thrust vector redirection due to roll angle 

perturbations, as shown in Figure IV. A. 1.-4. A comparison of Plots B. 1.-15, B. 1.-16 
shows the direct relationship of roll angle to the lateral acceleration. In the low speed region 
of the flight envelope, the yaw rate component, Y r ~ 0, due to the lack of forward velocity. 

This component tends to translate the forward velocity into the lateral velocity. As the 
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forward velocity increases this term directly increases in the form Y r ~ -V u . The geometry 
of the wing-tip RCS jets, as shown in Figure IV. A. 1.-4, relates the operations of the lateral 

stick to purely rolling motions. Assuming symmetry and the wing-tip RCS jet's lack of a 
direct effect on the lateral velocity, results in Y a - 0. The primary coupling of the yaw 

rotational motions to the lateral translations is provided by the tail-end RCS jets via rudder 
pedal control, Y^. Plot B. 1.-17 shows the direct relation of the rudder pedals to the lateral 

acceleration. This is due to the relatively large moment arm of the tail-end jet's physical 
configuration. 

The next group of simplifications are directed at the yaw/directional components. Using 
the geometry and symmetry arguments of above, the low speed characteristics of the 
couplings of the roll rate to the yaw rate is threw N p . The lateral velocity effects on the yaw 

rate are primarily due to the non-symmetric vehicle body configuration along the x axis (i.e. 
cross sectional area of the vertical stabilizer/rudder when compared to that of the 
forward/nose section). During low speed flight the yaw rate is not very sensitive to lateral 
translations and thus N v ~ 0. It is interesting to note that this term is important when dealing 

with a single main rotor helicopter with a tail rotor. This is due to the tail rotor's sensitivity 
to local sideslip and it's generally high main rotor disk loading [22]. 

The final area of simplification is directed at the roll dynamics. During low speed flight, 
the effects of small perturbations in lateral velocity on the roll rate, L v , can be neglected. 

This, however, is not the case for larger lateral velocities due to the size and orientation of 
the Harrier II AV-8B's critical wing. 
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Incorporating the above assumptions within the lateral-directional model of EQ(B.l.- 
18) results in the following low speed/powered-lift lateral model. 
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(B. 1.-19) 


An analysis will now be conducted to better understand the lateral responses of the low 
speed region of the flight envelope. This analysis will also provide some insight to further 
simplifications of the system model. The characteristic equation of the simplified system of 
EQ(B. 1.-19) is given by: 


a laT “ det 


S(S-Lp) -L r 0 
-SN p S - N r 0 
-Y* 0 S - Y v _ 


which results in 

A LAT = S(S - V (S ' N r )(S ‘ Y v) ' SL r N p< S ' Y v) 


(B. 1.-20) 


(B.l.-21a) 


or 


- S 4 - (L p +N r +Y v )S 3 + (L p N r +Y„L p +Y > N r -L,N p )S 2 + Y.^-L^S 


(B.l.-21b) 
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The classical lateral dynamics can be described by the characteristics equation: 


^LAT = (S + T s )(S + T r )(S 2 + 2w d 5 d S + w d 2 ) 


(B. 1.-22) 


where 


T s = The pole associated with the time constant of the spiral mode 

T r = The pole associated with the time constant of the roll subsidence mode 

w d = Natural frequency of the dutch role mode 

5 d = Damping ratio of the dutch role mode 

In general the time constant of the spiral mode is very long. For this reason and the fact that 
only the pilot frequencies are considered, the spiral mode can be reduced to an integrator 
model, thus: 

a'lat = S(S + T r )(S 2 + 2w d 5 d S + w d 2 ) (B.l.-23a) 

or 

a'lat = S 4 + S 3 (T r + 2w d 5 d ) + S 2 (2w d 5 d T r + w d 2 ) + ST r w d 2 (B. 1 .-23b) 

The terms of EQ(B.l.-23b) can be related to those of EQ(IV.a.l.b.-21b) as follows: 

-V N .- Y .= T . + 2w <A (B.l.-24a) 
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(B. 1 .-24b) 


VWV'WA - 2w <A T r + »d 2 

Y.tLjNp-LpN,) =T r w d 2 (B.1.-24C) 

Plot B. 1.-18 shows both the long and short period modes of the yaw rate response due to 
an impulse on the rudder pedals. This plot shows a divergent dutch roll mode that 
maintains a relatively long period. From this, one can assume: 

w d 2 is small (B. 1.-25) 


Including this assumption within EQ(B.l.-24c) we have: 

Y v ~ 0 (B.l.-26a) 

^Np-0 (B.l.-26b) 


The resulting low speed/powered-lift lateral dynamics model for pilot frequencies is given 
by: 
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(B. 1.-27) 


The characteristic equation of this system is given by: 
a"lat = S 4 - S 3 (Lp+N t ) + S 2 L p N, 


(B. 1.-28) 
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B.2. - Identification of Vehicle Model 


The high order state space models have been developed to facilitate direct parameter 
identification techniques. This section will derive the relationships that are used to identify 
the model parameters and finally obtain a parameter set from an analysis of the responses of 
the simulation program. The parameter identification techniques that are considered, rely on 
transfer function representations. To interface with these identification approaches, a set of 
transfer functions will be derived from the system models developed above. The transfer 
function representations are obtained by applying the Laplacian techniques used in the 
previous sections. The key to this approach is the development of a set of coupling 
numerators. These numerators will provide the transmission path characteristics from the 
cockpit control mechanisms to a specific state variable. The transfer function of a specific 
control function is obtained by introducing the appropriate characteristic equation as the 
denominator. To simplify this analysis, the longitudinal and lateral dynamics will again be 
decoupled to the forms of EQ(B. 1.-16) and EQ(B. 1.-27). 


The longitudinal transfer functions are based on the state system of EQ(B. 1.-16) and 

the characteristics equation of EQ(B.l.-17a). The parameters that require identification 
within the longitudinal model are: M e , M q , X 0 , Xj, Zj, and Z w . The coupling numerator of 

the longitudinal stick to the pitch angle is given by: 


6(S) ' 

S.(S) 


LONG 


(S) = det 
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0 o 

Me 0 
0 S-Z w 


M e S(S-Z w ) 


(B.2.-1) 


The resulting transfer function is: 
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(B.2.-2) 


0(S) _ Me S _ H 
*.(S) ' a L 0N0 (S) = scs - M ,> 


The pitch rate transfer function is therefore: 


PCS) _ Me 

s.(S) = 


(B.2.-3) 


This transfer function is a reasonably accurate model of the short period response 
characteristics of the pitch rate in Plot B.2.-1. The the parameters M e and M q are directly 

identifiable from EQ(B.2.-3) and Plot B.2.-1. 

M e ~ 0.7 1 (B.2.-4a) 

M q ~-1.56 (B.2.-4b) 


The coupling numerator of the nozzle angle to the forward velocity is given by: 


V»(S) 

8j(S) 


A LONG^) ~ ^et 
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-0 


o 0 

S(S-M q ) 0 

0 s-z^ 


XjS(S-M q )(S - Z w ) 


(B.2.-5) 


The resulting transfer function is given by: 

VS) _ X,S(S-M q )(S ■ Z w ) _ 3t (B 2.-6) 

8j(S) A lqng^^ 

This provides a fairly accurate approximation to the initial time response of Plot B.2.-2. 
The parameter X: is obtained from the step height and is given by: 
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Xj - -0.017 


(B.2.-7) 


The pitch angle coupling term, Xq, can be directly obtained from gravitational components. 

Due to the low speed assumption, the pitching component will tend to follow the response 
characteristics of the nozzle angle, which results in: 

X 0 = -0.017 (B.2.-8) 

This is a good assumption if fixed nozzle angles are considered. 


The coupling numerator for the throttle to the vertical rate is given by: 
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(B.2.-10) 


The resulting transfer function is given by: 

V,(S) Z^S-Mq) Zr 

Sr® WS) ' (S ‘ Z »> 

Plot B.2.-3 shows a damped sinusoidal response instead of the first order response. This is 
due to the engine response characteristics as shown in Plot B.2.-4. Comparing these plots, 
the terms Zj. and Z w can be extracted and result in: 

Zp- 0.024 (B.2.-12a) 

Z w ~ -0.83 (B.2.-12b) 
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The lateral transfer functions are based on the state system of EQ(B. 1.-27) and the 

characteristics equation of EQ(B. 1.-28). The parameters that require identification in the 
lateral-directional model are: L p , N , N r , Y p , L a , L r , N ff , and Y n . The coupling numerator of 

the lateral stick to the pitch angle is given by: 


!^-A LAT (S)=det 
5 a (S) LAT 
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(S-N r ) 

0 


0 

0 

s 


L a S(S-N r ) 


(B.2.-13) 


The resulting transfer function is given by: 

»(S) L a S(S-N f ) L a 

8,(S) 4 lat (S) ' S < S V 

The roll rate transfer function is therefore: 


(B. 2.-14) 



(B.2.-15) 


This transfer function is a reasonable approximation to the short period response of Plot 
B.2.-5. The parameters L # and L p can be directly obtained from Plot B.2.-5 and EQ(B.2.- 

15): 


L a ~ 1.89 (B.2.-16a) 

L ~ -3.45 (B.2.-16b) 

p 


The coupling numerator for the lateral stick to the yaw rate is given by: 
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^A, AT (S) = det 
5 a (S) LAT 


N a S 2 (S - L p ) 


(B. 2.-17) 
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The resulting transfer function is given by: 

r(S ) N a S ; (S - L p ) N, 
5.(S) A LAT (S) ( s - N ,> 


(B. 2.-18) 


The parameter N a can be obtained by comparing the short term response of Plot B.2.-6 to 
EQ(B.2.-18). 


N a ~ 0.2 (B.2.-19) 

The term N r is not easily derived from the relation of EQ(B.2.-18) ot Plot B.2.-7. Plot 
B.2.-7 will be utilized to determine N r in a later examination. 


The coupling numerator for the lateral stick to the lateral velocity is given by: 
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(B. 2.-21) 


The resulting transfer function is given by: 

V v (S) L a Y<, L a Y, 
5 a (S) a' lat (S) S 2 (S-Lp) 


(B. 2.-22) 
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The above transfer function provides a manner in which the pitch coupling term, Y^, can 

be obtained. The difficulty with this approach is the presence of the double integrator. A 
more straight forward approach is to directly obtain from a comparison of Plots B. 1 15 

and B. 1.-16 which yields: 

Y^ ~ 0.56 (B.2.-23) 


The coupling numerator of rudder pedals to the yaw rate is given by: 


^a’lat(S) = det 
8 r (S) LAT 


S(S - L p ) 



K 0 

Nrr 0 
S 


= N rT S 2 (S - 


(B.2.-24) 


The resulting transfer function is given by: 

r(S) N„S 2 (S - L p ) N„ 
5 r (S) A’ LAT (S) " 


(B.2.-25) 


A comparison of EQ(B.2.-25) and Plot B.2.-7 shows a much more direct access to the 
parameter N r . The parameters, N r and are given by: 

N r ~ -0.83 (B.2.-26a) 

N - 0.5 (B.2.-26b) 

IT 


The coupling numerator for the rudder pedals to the roll angle is given by: 
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= L n S(S-N t ) 


(B.2.-27) 
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(S-N r ) 

0 S J 


The resulting transfer function is given by: 

<D(S) _ Lrr S ( S - N r) _ K 

5 r (S) A laT (S) S(S-L p ) 

The roll rate transfer function is therefore: 


(B. 2.-28) 


2^1 = JS- (B.2.-29) 

5 r (S) 

The parameter L n can be obtained from Plot B.2.-8 and the use of the result of EQ(B.2.- 
29). 


L n ~-0.25 


(B. 2.-30) 


The full rank linearized state model used in the development of the OCM pilot can be 
seen in Figure B.2.-1. The identified parameters are listed below. 
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M e = 0.7 1 
M q = -1.56 
X } = -0.017 
X 0 = -0.017 
Zp = 0.024 
Z w = -0.83 

L a = 1.89 (B.2.-31) 

L p = -3.45 
L r = -0.25 
N r = -0.83 
N p = 0.21 
N„ = 0.5 
= 0.56 
Yfr = 00 
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APPENDIX C 

A USERS GUIDE TO THE OCM SOFTWARE 

This appendix serves as a user's guide for the OCM pilot software. An overview of the 
OCM software implementation is presented. The algorithm for configuring the OCM is 
illustrated. The procedure for operating the OCM within the VSRA environment is 
illustrated and examples are provided. 

C.l. - Overview of the OCM Software 

The OCM software is designed to be utilized as an active element of the NASA-VSRA 
simulation environment. The OCM software system consists of four software modules and 
a configuration file, as shown in the block diagram of Figure C. l.-l. The modules and file 
are defined as follows: 

1. OCM.LIST.NML - This file serves as the configuration file for the OCM. The 

continuous time internal reference model (F,G,H matrices), pilot delay (Td), noise 
parameters (Q u , R y ), system initial conditions, forcing function time constant, and OCM 
cost function weightings (Q,R,S) are contained in this file. The file is arranged in a 
"namelist" format. 

2. OCM_SETUP - This module initializes and configures the OCM environment. 
OCM initial condition, internal reference model, pilot data, and cost function weights are 
read from the OCM_LIST file. A discrete time representation of the internal reference 
model is generated. Optimal control gains are calculated by solving the steady state matrix 
Ricatti equation. Steady state Kalman gains are computed and covariance matrices are 
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initialized. The state transition matrix of the predictor is generated. All necessary 
information is loaded into specific COMMON regions for use by other OCM functions. 

3. OCM_TRAJ - This module generates the time based trajectory that the OCM_PILOT 
is to follow. The time referenced command sequence is integrated and bandlimited to 
provide a full rank command. The command sequence is independent of the sampling rate 
of the simulation environment. The integration computations require sample period 
information. 

4. OCM_PILOT - This module performs the active computations involved in the 
closed loop participation of the OCM within the VSRA simulation environment. These 
include measurements of the VSRA state, obtaining the command trajectory from the 
OCM_TRAJ module, noise model generation (via Box-Meuller approach), delay 
progression, measurement and update of the Kalman filter estimates, time advance 
predictions, and control input calculations. The control inputs arc then applied to the VSRA 
through the cockpit control mechanism variables. 

5. OCM_SUBS - This module contains a pool of utilities that simplify the organization 
and implementation of the other OCM modules. Some of the utility functions include: 
matrix and linear algebra operations, Kalman filter gain and covariance progression 
generators, optimal control solver and gain generator, and a continuous time to discrete 
time converter. 

To interface the OCM software to the VSRA simulation environment, the main VSRA 
driver program (VSRA_DRIVER) was modified to accommodate the OCM system. The 
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modifications involved that allocation of various common areas to support OCM 
operations, and the implantation of the OCM_SETUP and OCM_PILOT modules at 
specific points within the simulation initialization and primary execution loops. In addition, 
two new VSRA commands were introduced to handle the OCM initialization phases and the 
flight simulation operations involving the participating OCM. The use of these commands 
will be explained in greater detail in later sections. To accommodate the output of the OCM, 
the file writing code of the PLOTDATA.FOR subroutine was modified to include OCM 
variables within the unformatted output data file VSRA_POLY.PLT. 

C.2. - Installing the OCM 

The OCM software is contained on a VAX Files-1 1 formatted tape labeled "PITT". The 
following files must be recovered from the tape: 


1) OCM_SETUP.FOR 

1 1 )NAMELIST.NML 

2) OCM.PILOT.FOR 

12) OCM.LIST.VRT 

3) OCM.SUBS.FOR 

13) OCM_LIST.LAT 

4) OCM_TRAJ.FOR 

14) POOA.VRT 

5) PLOTDATA.FOR 

15) POOA.LAT 

6) VSRA_OCM.FOR 

16) OCM.PLT.VRT 

7) OCM.PLOT.FOR 

17) OCM_PLT.LAT 

8) V SRA_OCM.OPT 

18) prrri5i.DAT 

9) FTP.OLB 

19) prmio.DAT 

10) *.OBJ 



Files 1-6 correspond to the primary OCM software modules. File 7 is the modified plotting 
routine for the OCM. File 8 is the special linking configuration for the OCM. File 9 is the 
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VSRA library supplied by NASA-LEWIS. The file set *.OBJ (10) corresponds to the pool 
of pre-compiled VSRA modules used during the linking procedures by VSRA_OCM.OPT. 
File 1 1 is the VSRA configuration file. Files 12 and 13 are the OCM configuration files for 
vertical and lateral maneuvers, respectively. Files 14 and 15 are the VSRA setup command 
files for vertical and lateral maneuvers, respectively. Files 16 and 17 are the plotting 
configuration files for the vertical and lateral maneuvers, respectively. Files 18 and 19 are 
the unformatted data files of the pilot flight operations for the lateral and vertical 
maneuvers, respectively. 

Files 1-6 should be compiled with the VAX Debug function enabled (to be consistent 
with the VSRA format). The VMS command string for compiling is as follows: 

$ FORTRAN/DEBU G/NOOPT/CROS S_REF/CONT =99 filename.FOR 

where "filename.FOR" is the appropriate Fortran file from the above list. Linking 
operations are controlled by a modified version of VSRA_DRIVER.OPT 
(VSRA_OCM.OPT). The modifications incorporate the OCM software modules during the 
link process. The linking command string is given by: 

$ LINK/DEBUG VSRA_OCM/OPTJLIB:FTP/LIB 

where the device LIB: contains the library FTP.OLB. 

A previous version of the OCM software utilized three external libraries: 1) SLATECH, 
2) IMSL, 3) DISPLA. These libraries offer many routines (particularly the SLATECH 
routine RICSOL, that solves various versions of the matrix Ricatti equation) that simplify 
the generation of the optimal control gains, muscular system time constants, Kalman filter 
gains and covariance matrices. To comply with the requirements that the software generated 
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in this research be completely self-supporting, the IMSL and SLATECH libraries were 
removed and the OCM was fitted with comparable algorithms. The algorithms utilized by 
the present version of the OCM are based primarily on iterative/time-varying solutions and 
are therefore rather sluggish. The DISPLA library was, however, retained because the 
plotting packages supplied by NASA-LEWTS were supported by DISPLA. 

C.3. - Configuring the OCM 

The configuration of the OCM defines the vehicle under control, control objectives, and 
the pilot description. The principle operations involved in the configuration of the OCM 
pilot are summarized in the following: 

1. - Develop and insert the continuous time internal reference model of the vehicle 
under control into the OCM configuration file, OCM_LIST.NML. Appendix B illustrates 
the construction of the Harrier II AV-8B low speed/powered-lift, pilot frequency model that 
is supplied. Within the OCM_LIST.NML, the two dimensional arrays (FM,GM,HM) 
correspond to the state space representation matrices (A,B,C) or (F,G,H). Section C.7 of 
this appendix provides a listing of the primary variables used in the OCM. The model is 
dimensioned by the variables (NOCM,MOCM,LOCM) which correspond to the system 
order, number of inputs and outputs, respectively. The arrays, FM and GM, describe the 
vehicle dynamics while the measurement array, HM, is primarily dependent on the display 
configuration. The two dimensional array, WM, corresponds to the disturbance distribution 
matrix. The sampling period of the OCM execution (typically the sampling period of the 
simulation program) is selected and specified by the variable Tl. This variable is primarily 
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for use in OCM applications that require execution at rates other that the fundamental 
frequency of the simulation environment. 

2. - Determine and code the time based command sequence of the desired trajectory. 
The trajectory generating code resides in OCM_TRAJ. The version of OCM_TRAJ that has 
been supplied provides a simple implementation of the rate driven command sequencer and 
the rate integrating/bandlimiting full rank command generator. To use this strategy, a rate 
driven trajectory must be defined in the form of pulse trains on the appropriate rate 
commands. The pulse trains, in this case, are implemented by a sequence of "EF" 
statements creating a string of step functions that are overridden by the step occupying the 
present interval. With the trajectory specified, the time constant of the forcing function 
bandlimiting filter is selected and specified by the variable AFORCE in the configuration 
file OCM_LIST.NML. This time constant typically ranges between 1 .5 and 4.5 seconds 
[2]. A value of 2,5 seconds has been utilized in the implementation supplied. The user may 
wish to insert his own trajectory defining code or route the trajectory information to the 
OCM from some external process via this routine. It is important to note that modifications 
to the OCM_TRAJ routine will require that the OCM_TRAJ.FOR file be re-compiled and 
the total software system be re-linked. 

3. - The control objectives of the pilot's task is defined in the form of cost function 
weightings. The values utilized reflect the manner in which the pilot will respond to the 
given situation. The values will typically vary from task-to-task and from vehicle-to- 
vehicle. The cost function weights are loaded into the two dimensional arrays QOPT, 
ROPT, GOPT within OCM.LIST.NML. The array QOPT corresponds to the definition of 
the acceptable maximum deviations of the errors in vehicle attitude and orientation from that 
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of the trajectory. ROPT defines the maximum deflection of the cockpit control mechanisms 
and is usually a function of the vehicle (see section V.B.3.). The array GOPT is adjusted to 
obtain the desired muscular system time constants. This typically requires a degree of 
iterative adjustment. The values supplied in this version reflect a subjective analysis of the 
piloted flight data and a limited knowledge of the target configuration and control 
objectives. 

4. - The final step is the selection of the pilot's inherent parameters. The discrete nature 
of the OCM and the VSRA simulation environment requires that the pilot delay be 
implemented as a chain of sample delay periods arranged in a FIFO buffer. The length of 
the buffer is determined by the number of sample periods needed to achieve the delay. The 
variable NDEL of the OCM_LIST.NML file, specifies the number of sample 
periods/elements of the buffer, that the pilot delay occupies. The remnant model noise 
sources, observation and motor, are selected according to the control objectives, vehicle 
under control, and the display configurations. The arrays STVU and STVY correspond to 
the noise model variances, Q u and R y , of the motor and observation noises respectively. 

C.4. - Executing the OCM within the VSRA Simulation Environment 

Executing the OCM software within the VSRA is relatively simple. The OCM 
operations are broken into two separate functions: 1) Initialization and preparation of the 
OCM environment, 2) Execution of the VSRA with the OCM actively participating in the 
flight control loops. These operations are provided by two VSRA commands: 
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/SOCM - This command executes the software module OCM_SETUP. The OCM 
environment is configured and loaded into specific common regions. This 
operation is typically performed after the vehicle has been trimmed with the 
TRIM command. The user may wish to modify the OCM_SETUP routine 
to create an external file of the pilot configurations instead of loading the 
common regions. The user will, however, have to provide the necessary file 
reading and common region loading facilities (possibly by an additional 
VSRA command). Upon the completion of this command, the pilot's 
parameters are displayed to the user. Again, the user may wish to modify 
the OCM_SETUP to have the pilot's parameters loaded into an external file. 

/ROCM - This command executes the primary simulation of the VSRA and enables 
the OCM operations. This command is tailored after the DYNC command 
with the exception that is utilizes the cockpit control mechanism deflections 
of the OCM instead of the dynamic check tests. This command can only be 
executed after the use of the TRIM and SOCM commands. It is important to 
note that the OCM is designed to operate with ONLY the AV-8B aircraft 
dynamics and not those of the YAV-8B. 

As mentioned previously, the algorithms utilized by the present version of the OCM 
software are based on iterative/time-varying solutions of the matrix Ricatti equation [27,28] 
(instead of the SLATECH RICSOL routine). This causes some complications in the 
techniques utilized to generate the optimal control gains, muscular system time constants. 
These are typically in the form of trial-and-error iterations of the OCM_SETUP operations 
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by using the command SOCM until the desired muscular system time constants are 
obtained. 

The overall operation of the OCM within the VSRA can be summarized as follows: 

1. - Configure the OCM environment by preparing the OCM_LIST file according to 
procedures of section C.3 of this appendix. The VSRA simulation environment is 
configured by preparing the NAMELIST.NML file. 

2. - Enter the VSRA simulation environment and issue the command SOCM. This will 
generate the OCM control and estimation gains and loads the specific common areas. The 
values of the pilot parameters are displayed upon completion. The user should examine the 
pilot parameters to determine if a satisfactory pilot profile has been obtained. If so, proceed 
to step 3, if not, exit the VSRA and modify the OCM_LIST file, then repeat step 2. 

3. - With the proper pilot parameters resident within the VSRA environment, the flight 
simulation may begin. The user issues the command ROCM and the VSRA proceeds to 
execute the simulation according to the NAMELIST.NML file. The output of the VSRA is 
deposited in the unformatted data file VSRA_POLY.PLT. 

C.5. - Output Generated by the OCM 

The time based output sequence of the VSRA simulation is deposited in the unformatted 
data file VSRA_POLY.PLY. This file is generated by a modified version of the subroutine 
PLOTDATA.FOR. The modifications were made to accommodate the output of the OCM. 
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In addition to the OCM software that has been provided, a pair of unformatted data files 
containing the human piloted flight data (PITT151.DAT [lateral tracking maneuver] and 
PITT210.DAT [vertical tracking maneuver]) are also included. 

The data of the unformatted files can be plotted with the routine OCM_PLT.FOR. This 
routine is a modified version of the plotting package supplied by NASA-LEWIS and is 
based on the DISPLA library. This routine relies on the configuration file OCM_PLT.SRC 
to provide the necessary default plotting information. OCM_PLT.FOR permits the user to: 
1) plot the OCM output, 2) plot the piloted flight data (either lateral or vertical maneuvers), 
or 3) plot a comparison of both (as shown in OCM comparison plots). This routine can 
also be utilized to plot data from strictly DYNC runs. 


The plotting activities are arranged according to the cockpit control mechanisms and the 
dominant aircraft responses of those controls. These break-downs are given below: 


1) Longitudinal Stick 

Pitch angle 

Pitch rate 

Longitudinal stick 

2) Lateral Stick 

Roll angle 

Roll rate 

Lateral stick 


Lateral position Lateral velocity 

3) Rudder Pedals 

Yaw angle 

Yaw rate 

Rudder pedals 


Sideslip 
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4) Throttle 


Altitude 


Vertical rate Throttle 


Engine speed 

5) Nozzle Angle Forward velocity Nozzle angle 

Forward position 

To execute the OCM plotting routine issue the following VMS commands: 

$ FOR OCM.PLT.FOR 
$ LINK OCM_PLT,LIB :DISPLA/LIB 
$ RUN OCM.PLT 

The term LIB :D IS PL A/LIB corresponds to the link search of the DISPLA library residing 
on device LIB:. This command will depend on the VMS configuration and File structure 
being used. The routine will respond with the following question: 

PLOTS?? (Flight:0, OCM:l, Both:2) > 

The user should select the appropriate data set to plot and give a numerical answer. If the 
user selects either 0 or 2, the routine will request the desired flight maneuver. 

Task??(Vert:0, Lat:l)> 
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The routine will then request the control mechanism group to be plotted by: 

What kind of control input? 

1 = Longitudinal Stick 

2 = Lateral Stick 

3 = Rudder Pedals 

4 = Throttle 

5 = Nozzle 

Enter the appropriate number : > 

The routine will then read in the VSRA_POLY.PLT, PITT151.DAT, or PiiT210.DAT 
files and proceed to generate the appropriate plots. The plotted output of this routine will be 
a set of plot files that are in the Tektronix 4010 graphics format. 

C.6. - Examples of VSRA/OCM Execution 

As a conclusion to this overview of the OCM software, an example will now be 
presented. This example illustrates the configuration and execution of the OCM within the 
VSRA environment. The listings of the this interactive session was captured via terminal 
monitoring facilities. The maneuver in question is the vertical tracking maneuver. 

Vertical Tracking Maneuver Example 

The OCM.LIST.NML file for this maneuver is provided in the file OCM_LIST.VRT. 
The VSRA setup file is P00A.VRT. Both of these files must be copied to their operational 
names of OCM_LIST.NML and POOA.COM, respectively. 
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$ COPY OCM_LIST. VRT OCM_LIST.NML 

Copying USR5: [02350. MZVSTOL. HARRIER. SPOOL] OCM_LIST. VRT;6 to USR5: [02350. MZVSTOL 
. HARRIER. SPOOL] OCM_LIST. NHL; 131 3 blocks 

$ COPY P00A.VRT POOA.COM 

Copying USR5: [02350. MZVSTOL. HARRIER. SPOOL] POOA. VRT; 3 to USR5: [02350. MZVSTOL. HAR 
RIER. SPOOL] POOA.COM; 25 2 blocks 

$ RUN VSRA_OCM 

VAX DEBUG Version V5 . 0~00 MP 

DEBUG- I -INI UAL, language is FORTRAN, module set to VSRA_OCM 
DBG> @P00A 


WELCOME TO THE VSRA VAX SIMULATION PROGRAM 


Which aerodynamics would you like to use ? 

1 ■ AV-8B aerodynamics 

2 ■ YAV-8B aerodynamics 

Enter the appropriate number: 1 


ENTER ?? FOR COMMAND LIST 


VSRA> ?? 


VSRA SIMULATION INTERACTIVE COMMAND LIST: 


/TRIM RUN TRIM PROGRAM 

/TRIM XX RUN TRIM PROGRAM FOR XX CYCLES 

/PRNT X TO PRINT COMMON BLOCKS: X-2 FOR XFLOAT COMMON BLOCK 

X-3 FOR USER COMMON BLOCKS 
X-A FOR IFIXED COMMON BLOCK 

/DATA INPUT DATA TO XFLOAT AND I FLOAT FROM DATA FILE 

/UDAT INPUT DATA TO USER COMMON BLOCKS FROM DATA FILE 

/I CRN XX EXECUTE I.C. RUN FOR XX CYCLES 

/DYNC X RUN DYNAMIC CHECK - X>0 PRINT HEADER AND DATA FOR X SECS 

/OPRN XX TO EXECUTE RUN FOR XX SECONDS 

/MESS SEND MESSAGE TO PRINT OUTPUT FILE 

/CHNG TO CHANGE COMMON BLOCK VARIABLE INTERACTIVELY 

/STAB TO COMPUTE STABILITY DERIVATIVES 

/FGTB TO WRITE F & G MATRICES 

/SAVT TO SAVE TRIM VALUES 

/REST TO RESTORE TRIM VALUES 

/CHGC TO CHANGE CONFIGURATION ID 

/SOCM TO INITIALIZE AND SETUP THE OCM 

/ROCM TO EXECUTE THE OCM WITHIN THE VSRA 

/END TO TERMINATE THE PROGRAM 

?? TO PRINT THIS LIST 
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VSRA> /SOCM 


****** * ******** * *********** * * * 
* * 

* OCM SETUP * 

* * 
****************************** 


DISCRETE TIME REPRESENTATION COMPLETED 
OPTIMAL CONTROL GENERATED 
KALMAN FILTER GENERATED 
TIME ADVANCE PREDICTOR GENERATED 


OCM SETUP COMPLETED 


PILOT PARAMETERS 
PILOT DELAY 0.20 

FORCING FUNCTION TIME CONSTANT 2.500 


MUSCULAR SYSTEM TIME CONSTANTS 


CONTROL 

MECHANISM 

i 

0.15929 

CONTROL 

MECHANISM 

2 

0.00000 

CONTROL 

MECHANISM 

3 

0.15389 

CONTROL 

MECHANISM 

A 

0.15397 

CONTROL 

MECHANISM 

5 

0.13A71 

OR NOISE VARIANCES 



CONTROL 

MECHANISM 

1 

1.00000 

CONTROL 

MECHANISM 

2 

1 . 00000 

CONTROL 

MECHANISM 

3 

20.00000 

CONTROL 

MECHANISM 

A 

0.50000 

CONTROL 

MECHANISM 

5 

0.50000 


OBSERVATION NOISE VARIANCES 


DISPLAYED 

VARIABLE 

i 

0.85000 

DISPLAYED 

VARIABLE 

2 

0.60000 

DISPLAYED 

VARIABLE 

3 

1.00000 

DISPLAYED 

VARIABLE 

A 

0.71000 

DISPLAYED 

VARIABLE 

5 

0.85000 

DISPLAYED 

VARIABLE 

6 

0.60000 

DISPLAYED 

VARIABLE 

7 

1.00000 

DISPLAYED 

VARIABLE 

8 

0.71000 

DISPLAYED 

VARIABLE 

9 

1.00000 

DISPLAYED 

VARIABLE 

10 

0.71000 

DISPLAYED 

VARIABLE 

11 

0.85000 

DISPLAYED 

VARIABLE 

12 

0.60000 



VSRA> /TRIM 


****************************** 
* * 

* TRIM MODE * 

* * 

* * y? * ★ vt * ^ ★ yc ^ yt * yr yr ^ vt vt y^r yir yir vt y*c 


TRIM IS SUCCESSFUL AFTER 205 CYCLES 
VSRA> /ROCM 


yt yr yc * y« ★ * * * fr/c * * * y? * * * * Vr ** y^r * yt * * * * * 

* * 

* OCM EXECUTION * 

y k * 

yc * * ******* * **** * * * ** * ******** * 

VSRA> /END 

VSR A SIMULATION COMPLETED SUCCESSFULLY 

%DEBUG-I-EXITSTATUS, is ' "4SYSTEM-S-N0RMAL, normal successful completion' 
DBG> EXIT 
$ DIR * . PLT 

Directory USR5 : [02350. MZVSTOL. HARRIER. SPOOL] 

SMAR THIS. PLT; 1 TRIMHIS.PLT; 1 VSRA_POLY.PLT; 1 

Total of 3 files. 

$ 
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C.7. - Listing of the OCM variables 


The OCM data structure has been implemented with the user in mind. No high 
performance array declaration or common area structures have been constructed. This will 
allow the user to best tailor the OCM environment to his application. The following list 
defines the primary variables utilized within the OCM. 

£ 

- Dimension of the OCM's internal reference model of the system under 
control. 

- Number of cockpit control mechanisms (inputs). 

- Number of displayed variables (system displays). 

- Number of simulation sample periods per pilot delay 

REALIS 

T1 - Sampling period of the OCM execution in seconds. 

X(*) - System state vector. 

U(*) - Vector of cockpit control mechanisms 

FM(*,*) - System matrix of the internal reference model. 

GM(*,*) - Input distribution matrix of the internal reference model. 

HM(*,*) - Measurement matrix of the internal reference model 

WM(*,*) - Disturbance distribution matrix 

PHIM(*,*) - Discrete time system matrix 

GMA(*,*) - Discrete time input distribution matrix 

U0(*) - Initial conditions of the cockpit control mechanisms. 

UC(*) - Desired control vector. 

UCKDZ(*,*) - FIFO buffer for the control input delay. 
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YKDZ(*,*) 

QOPT(*,*) 

ROPT(*,*) 

GOPT(*,*) 

KOPT(*,*) 

POPT(V) 

PMOPT(*,*) 

PHI0(*,*) 

GMA0(*,*) 

Q10PT(*,*) 

STVU(*) 

STVY(*) 

VUM(*,*) 

VYM(*,*) 

PKAL(*,*) 

PMKAL(*,*) 

PHIK10(* ,*) 

GMAK00(V 

GMAK10(V 

Kl(*,*) 

H1(V) 

PHI1(*,*) 

GMA1(*,*) 


- FIFO buffer for the system output delay. 

- Cost function weights for the system states 

- Cost function weights for the control inputs 

- Cost function weights for the control rates. 

- Optimal control gain state feedback matrix 

- Update matrix of the optimal control gain generator 

- Measurement update matrix of the optimal control gain generator 

- Augmented system matrix for the solution of the optimal control gains. 

- Augmented input distribution matrix for the solution of the optimal control 
gains. 

- Augmented cost function state weighting matrix for the solution of the 
optimal control gains. 

- Variances of the motor noise sources. 

- Variances of the observation noise sources. 

- Covariance matrix of the motor noises. 

- Covariance matrix of the motor noises. 

- Covariance matrix of the time-varying Kalman filter solution 

- Measurement covariance matrix of the time-varying Kalman filter solution 

- Full rank optimal control gain matrix 

) - State feedback optimal control gains 

) - Muscular system optimal control gains (time constants) 

- Kalman filter output error correction gain matrix 

- Augmented measurement matrix for the solution of the Kalman filter 
gains. 

- Augmented system matrix for the solution of the Kalman filter gains. 

- Augmented input distribution matrix for the solution of the Kalman filter 
gains. 
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Q1(V) 

Xl(*) 

X1MINUS(*) 

FM1(*,*) 

K2(*,*) 

ZKDZ(* ,*) 

X2(*) 

ZETA(*) 

XR(*) 

XRKD(*) 

XR1(*) 

AR1(* ,*) 

BR1(*,*) 

AFORCE 

XR2(*,*) 

AR2(*,*) 

BR2(*,*) 

XR3(*,*) 

IDNT(* *) 

NULL(*,*) 

W(*,*) 

Wl(* ,*) 
W2(*,*) 
WINV(*,*) 


- Augmented noise model covariance matrix for the solution of the Kalman 
filter gains. 

- Augmented estimate of the state vector from the Kalman filter. 

- Update estimate of the augmented state vector of the Kalman filter. 

- Augmented continuous time system matrix for the generation of the state 
transition matrix of the least-squared predictor. 

- State transition matrix of the least-squared predictor. 

- FIFO buffer for the state estimate delay. 

- Time advanced prediction of the system augmented state 

- Predicted state estimate 

- Trajectory reference state vector. 

- Delayed trajectory reference state vector 

- State vector of the bandlimiting filters of the trajectory generators 

- System matrix of the bandlimiting filters 

- Input distribution matrix of the bandlimiting filters 

- Time constant of the forcing function bandlimiting filter 

- State vector of the integrated and distributed command sequence 

- System matrix of the command sequence integrator 

- Input distribution matrix of the command sequence integrator 

- State vector of the rate command sequence 

- Utility identity matrix 

- Utility zero matrix. 

- Working array 

- Working array 

- Working array 

- Working inversion array 
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Figure LA.-l - Block Diagram of the closed loop piloted control structure 


COMMAND 


COCKPIT CONTROL 
MECHANISMS 


ATTITUDES 

AND 

ORIENTATIONS 



Figure II.- 1. - Block diagram of the basic human controller characteristics 
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Figure II.A.-l. - Block diagram of the internal structure of the McRuer-Krendel model 
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Figure II.B.-l. - Simple block diagram of the OCM within a control environment. 


125 













secondary 

RESPONSES 


PRIMARY 

RESPONSE 


COCKPIT CONTROLS AIRCRAFT DYNAMICS VEHICLE RESPONSES 


Figure IIT.B.l.-l. - A block diagram of the separation of the primary and secondary 
response characteristics associated with the decoupled transfer function models. 
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Figure III.B.2.-1. - Control loop closing strategy for the design of single variable pilot 

mechanisms 
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Figure III.B.2.-2. - Illustration of the destabilizing distortions associated with a delay 

Root-Locus 



Figure in. B. 2.-3. - Illustration of a delay approximation using a large pole set 


127 


ORIGINAL PAGE IS 
OF POOR QUALITY 


comwo 

SET 



Figure III.B.3.-1. - Multi-variable control structure for the insertion of the McRuer- 

Krendel pilot models 



Figure III.B.3.-2. - Multi-variable control structure for executing simple flight control 

maneuvers. 
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Figure IV.A.l.-l. - Illustration of the modification of the engine thrust vector due to a 
reduction in nozzle angle during low speed/powered-lift activities 



Figure IV. A. 1.-2. - Illustration of the modification of the engine thrust vector due to an 
increase in engine speed during low speed/powered-lift activities 
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FORWARD RCS JET 
EXHA l’ ST VECTOR 



Figure IV. A. 1.-3 - Illustration of the rotational motion due to the operation of the 
forward RCS jet from the backward deflection of the longitudinal stick during low 

speed/powered- lift activities. 
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WING-TIP RCS JET 
EXHAUST VECTOR 



Figure IV. A. 1.-4. - Illustration of the rotational motion due to the operation of the 
wing-tip RCS jets from the deflection of the lateral stick during low speed/powered-lift 

activities 
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Figure IV. A. 1.-5. - Illustration of the rotational motion due to the operation of the tail- 
end RCS jet from the deflection of the rudder pedals during low speed/powered-lift 

activities 



Figure V.A.l.-l. - Root-Locus of the pitch control pilot mechanism 
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Figure V.B.l.-l. - Illustration of the target orientation and vehicle motion during the 
Vertical Tracking Hover maneuver 



Figure V.B.1.-2. - Illustration of the target orientation and vehicle motion during the 

Lateral Tracking Hover maneuver 
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Figure V.B.4.-1. - Diagram of a trajectory reference generator for the pitch angle 
components driven by a pitch rate command sequence. 
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Figure V.B.4.-2. - Vertical rate and yaw rate command sequence for the vertical 

tracking maneuver. 
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Figure A.l. -1. - Block diagram of the internal structure of the OCM 
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Figure C.l. - 1. - Block diagram of the software modules and configuration files of the 

OCM simulation environment. 
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TABLE LISTING 


xx u 


— Longitudinal position- 

V u 


Forward velocity 

0 


Pitch angle 

q 


Pitch rate 

ZZw 


Altitude 
Vertical rate 

V 



w 

9 


Roll angle 
Roll rate 

p 



V 


Yaw angle 
Yaw rate 

r 



YY V 


Lateral position 
- Lateral velocity -* 

. v v _ 



Table IV.A.3.-1 - List of the state vector variables for the high order Harrier model. 
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-5 e - 

5 t 


5 


a 



- Longitudinal stick position (inches) - 
Nozzle angle control position (degrees) 

Power level/throttle position (percent) 
Lateral stick position (inches) 

- Rudder pedal position (inches) - 1 


Table IV.A.3.-2. - List of the control vector variables for the high order Harrier model. 


V 

(knots) 

w e 

(rads/sec) 

— 

5e 

k to 

(in-sec 3 )' 1 

K x 

(in-sec 2 ) 1 

20 

4.7 

0.28 

0.041 

2.32 

| 40 

4.9 

0.19 

0.045 

2.59 

i 

60 

5.1 

0.14 

0.051 

2.91 

80 

5.4 

0.10 

0.049 

2.83 

! 100 

S i 

5.7 

0.07 

0.048 

2.77 


Table IV.B.L-1. - Transfer function parameters for the pitch response of the Black 

Hawk. 
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V 

k alt 

I 

a ALT 1 

(knots) 

(. fcet 2 ) 
in-sec z 

(sec)' 1 ' 




20 

6.6 

0.30 | 

j 

! 40 ; 

’ l 

7.0 

| 

0.45 ; 

» 

1 60 

1 

7.0 

0.50 

i 

80 

7.3 

0.50 

j 100 

i J 

7.5 

i 

0.50 i 

j 


Table TV. B. 1.-2. - Transfer function parameters for the altitude response components 

of the Black Hawk. 
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Altitude Control 

| 

Vertical Rate Control 
I Velocity Control 

t 

l 

Roll Control 
Lateral Velocity Control 

j 

Heading Control 
; Yaw Rate Control 
Sideslip Regulation 


52 - 


36 - & 



-14. ( 


deg 


knot-sec 

20 ( — ) 

‘ Meg-sec' 

in 


*1 


0.0 0.0 


3.6 8.0 


l - 1 (^) j 0.8 6.0 ! 


26. (x- 12 — ■) 0.8 

Meg-sec 7 




32- (a^=r) 0.8 20. j 


Meg-sec 


20 . 


0.8 0.0 1 


J 


Table V.A.l.-l. - Parameter list of the additional Harrier pilot mechanisms. 
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V 

(knots) 

l 1 

( in ) 
Meg-sec 7 





T 


20 : 

1.31 

\ 

40 | 

1.48 

60 

1.53 

80 j 

2.07 

100 

1.82 


Table V. A. 1.-2. - Parameters of the pitch control pilot of the Black Hawk 
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A = 0.8 (sec) 1 



i 


20 

; 5.77 

5 

40 

5.68 

60 

! 4.97 

t 

80 

| 5.07 

t 

j 

100 

! 5.19 


B = 5.5 (sec) 1 


Table V.A.1.-3. - Parameters of the altitude control pilot for the Black Hawk 



A = 0.8 (sec)' 1 


K 

in ^ 

© B = 0.0 (sec)* 1 


20 

1.82 

40 

1.74 

60 

1.43 

80 

1.51 

100 

j 1.57 
i 


Table V. A. 1.-4. - Parameters of the altitude rate control pilot for the Black Hawk 


V 

(knots) 
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V 

K 

(knots) 

( in ) 

Meg- sec' 



20 

! 

0.57 | 

j 

40 

0.57 1 

i 

60 ! 

i 

0.68 

80 

0.68 

100 

0.68 

j 


A = 0.8 (sec) 1 
B = 0.0 (sec) 1 


Table V. A. 1.-5. - Parameters of the roll angle control pilot for the Black Hawk 



V 

K 

A 

B 

(knots) 

( in ) 1 
Meg-sec 7 ; 

(sec) 1 

(sec) 4 





20 

0.63 

2.0 

i 

1.5 i 

• 40 

i 

o.6o ; 

l- 4 ! 

| 

0.75 . 

! 1 

; 60 

i 

0.43 

0.8 

o 

o 

80 

0.51 

i 

0, 

0.0 

; 100 

0.65 

0.8 

1 

i 0.0 

i . . 


Table V.A. 1 .-6. - Parameters of the heading control pilot for the Black Hawk 
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V 

(knots) 

K ! 
111 ) 

v deg-sec / 

I i . . . 

— — J 




I 


\ 



I 20 

-0.25 

— 

; 40 

-0.45 

— 

! 60 

1 

-0.55 


o 

oo 

1 -0.55 

! 


1 

j 100 

j 

\ -0.55 

i 

I 

1 


A = 0.8 (sec)' 1 
B = 0.0 (sec)' 1 


Table V. A. 1.-7. - Parameters of the Sideslip regulation pilot of the Black Hawk 
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cockpit 

control 

^"vmechanism 

maneuver 

Longitudinal 

Stick 

Lateral 

Stick 

Rudder 

Pedals 

Throttle 

Nozzle 

Angle 

Pitch 

Reorientation 

Pitch 

Control 

Pilot 

Roll 

Control 

Pilot 

Heading 

Control 

Pilot 

Altitude 

Control 

Pilot 

Velocity 

Control 

Pilot 

Velocity 

Translation 

Pitch 

Control 

Pilot 

Roll 

Control 

Pilot 

Heading 

Control 

Pilot 

Altitude 

Control 

Pilot 

Velocity 

Control 

Pilot 

Altitude 

Translation 

Pitch 

Control 

Pilot 

Roll 

Control 

Pilot 

Heading 

Control 

Pilot 

Altitude 

Control 

Pilot 

Velocity 

Control 

Pilot 

Heading 

Modification 

Pitch 

Control 

Pilot 

Roll 

Control 

Pilot 

Heading 

Control 

Pilot 

Altitude 

Control 

Pilot 

Velocity 

Control 

Pilot 

Altitude Rate 
Translation 

Pitch 

Control 

Pilot 

Roll 

Control 

Pilot 

Heading 

Control 

Pilot 

Altitude 

Rate 

Pilot 

Velocity 

Control 

Pilot 

Flat Turn 

Pitch 

Control 

Pilot 

Roll 

Control 

Pilot 

Yaw Rate 

Control 

Pilot 

Altitude 

Control 

Pilot 

Velocity 

Control 

Pilot 

Coordinated 

Turn 

Pitch 

Control 

Pilot 

Roll 

Control 

Pilot 

Sideslip 

Regulation 

Pilot 

Altitude 

Control 

Pilot 

Velocity 

Control 

Pilot 


Table V. A. 2.-1. - Table of various flight control maneuvers 

and their associated configurations of 
Harrier SVPMs. 
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cockpit 

control 

^..mechanism 

maneuver 

Longitudinal 

Cyclic 

Stick 

Lateral 

Cyclic 

Stick 

Rudder 

Pedals 

Collective 

Stick 

Pitch 

Reorientation 

Pitch 

Control 

Pilot 

Roll 

Control 

Pilot 

Heading 

Control 

Pilot 

Altitude 

Control 

Pilot 

Altitude 

Translation 

Pitch 

Control 

Pilot 

Roll 

Control 

Pilot 

Heading 

Control 

Pilot 

Altitude 

Control 

Pilot 

Heading 

Modification 

Pitch 

Control 

Pilot 

Roll 

Control 

Pilot 

Heading 

Control 

Pilot 

Altitude 

Control 

Pilot 

Altitude Rate 
Translation 

Pitch 

Control 

Pilot 

Roll 

Control 

Pilot 

Heading 

Control 

Pilot 

Altitude 

Rate 

Pilot 

Coordinated 

Turn 

Pitch 

Control 

Pilot 

Roll 

Control 

Pilot 

Sideslip 

Regulation 

Pilot 

Altitude 

Control 

Pilot 


Table V.A.2.-2. - Table of various flight control maneuvers 

and their associated configurations of 
Black Hawk SVPMs. 




VARIABLE 

TRIMMED 

COMMAND 

Pitch (9 C ) 

7.3 degrees 

17.3 degrees 

Roll (<j> c ) 

0.0 degrees 

0.0 degrees 

i 

Heading (\|/ c ) 

0.0 degrees 

0.0 degrees 

Altitude 

100 feet 

100 feet 

Velocity 

10 knots 

; 10 knots 

i 


Table V.A.3.-1. - Pitch reorientation command sequence for the multi-variable Harrier 

pilot. 
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VARIABLE 

TRIMMED 

COMMAND 

i 

[ 

! Pitch (9 C ) 

' 

6.5 degrees 

6.5 degrees 

Roll (4> c ) 

0.0 degrees 

0.0 degrees 

Heading (\)/ c ) 

i 

1 

- 

0.0 degrees 

| 0.0 degrees 

i 

i 

Altitude 

100 feet 

100 feet ; 

Velocity 

• J 

1 25 knots 

j 20 knots 


Table V.A.3.-2. - Velocity translation command sequence for the multi-variable Harrier 

pilot. 
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— 

j VARIABLE 

J 

TRIMMED ! 

1 

! 

J 

COMMAND 


Pitch (0 C ) 

\ 

\ 

2.8 degrees j 

2.8 degrees 

t 

— - 

Roll «t> c ) 

1 0.0 degrees 

20.0 degrees 

T- 

Sideslip (P c ) 

i 

| 0.0 degrees 

i 

0.0 degrees > 

i 

i 


: Altitude 

; 200 feet 

| j 

t 

200 feet 


Table V.A.3.-3. - Coordinated turn command sequence for the multi-variable Black 

Hawk pilot 
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Table V.B.3.-2. - Cost function weights of the OCM for the lateral tracking maneuver 
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Table V.B.3.-3. - Magnitudes of the OCM motor noise sources applied to each cockpit 

control mechanism 
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Table V.B.4.-1. - Trimmed values of the Harrier simulation environment for the 
vertical tracking precision hover maneuver. 
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Table V.B.4.-2. - Trimmed values of the Harrier simulation environment for the lateral 

tracking precision hover maneuver. 
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PLOT GROUP 
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TIME (SEC) 

Plot IV.A.2.-1. - Pitch rate response of the Harrier AV-8B due to a unit impulse 
deflection of the longitudinal stick in a near hover 



TIME (SEC) 


Plot IV. A. 2.-2. - Roll rate response of the Harrier AV-8B due to a unit impulse 
deflection of the lateral stick in a near hover 
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TIME (SEC) 

Plot IV. A. 2.-3. - Forward velocity response of the Harrier AV-8B due to a 5 degree 
impulse deflection of the nozzle angle at 10 knots 


TIME (SEC) 

Plot IV.A.3.-1. - Long term pitch rate response of the Harrier AV-8B due to a unit 
impulse deflection of the longitudinal stick at 10 knots 
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Plot IV. B. I.-l. - Pitch rate response of the Black Hawk UH-60A due to an impulse of 
the longitudinal cyclic stick at 60 knots. 
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Plot IV.B.1.-2. - Pitch rate response of the Black Hawk UH-60A due to an impulse of 
the longitudinal cyclic stick at 100 knots. 
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Plot IV.B.1.-3. - Roll rate frequency response of the Black Hawk UH-60A due to the 
operation of the lateral cyclic stick at 80 knots. 
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Plot IV. B. 1.-4. - Altitude rate response of the Black Hawk UH-60A due to a step of 
the main rotor collective stick at 60 knots. 
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Plot V.A.3.-1. - Pitch angle response during a pitch reorientation maneuver by a multi- 
variable McRuer-Krendel pilot inserted in the Harrier AV-8B 



Plot V.A.3.-2. - Longitudinal stick response of a multi-variable McRuer-Krendel pilot 
during a pitch reorientation maneuver in the Harrier AV-8B 
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Plot V. A. 3.-3. - Forward velocity response during a pitch reorientation maneuver by a 
multi-variable McRuer-faendel pilot inserted in the Harrier AV-8B 



Plot V.A.3.-4. - Nozzle angle response of a multi-variable McRuer-Krendel pilot 
during a pitch reorientation maneuver in the Harrier AV-8B 
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Plot V.A.3.-5. - Altitude response during a pitch reorientation maneuver by a multi 
variable McRuer-Krendel pilot inserted in the Harrier AV-8B 
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Plot V.A.3.-6. - Throttle response of a multi-variable McRuer-Krendel pilot during a 
pitch reorientation maneuver in the Harrier AV-8B 
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Plot V.A.3.-7. - Forward velocity response during a velocity translation maneuver by a 
multi-variable McRuer-Krendel pilot inserted in the Harrier AV-8B 



Plot V.A.3.-8. - Nozzle angle response of a multi-variable McRuer-Krendel pilot 
during a velocity translation maneuver in the Harrier AV-8B 
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Plot V.A.3.-9. - Pitch angle response during a velocity translation maneuver by a multi- 
variable McRuer-Krendel pilot inserted in the Harrier AV-8B 
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Plot V. A. 3.-1 1. - Altitude response during a velocity translation maneuver by a multi- 
variable McRuer-Krendel pilot inserted in the Harrier AV-8B 
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Plot V.A.3.-13. - Roll angle response during a coordinated turn maneuver by a multi- 
variable McRuer-Krendel pilot inserted in the Black Hawk UH-60A 
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Plot V.A.3.-14. - Yaw angle response during a coordinated turn maneuver by a multi- 
variable McRuer-Krendel pilot inserted in the Black Hawk UH-60A 
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Plot V.A.3.-15. - Sideslip angle response during a coordinated tum maneuver by a 
multi-variable McRuer-Krendel pilot inserted in the Black Hawk UH-60A 


ALTITUDE POSITION RESPONSE 



^UT kFEET) VS (S£C> 


Plot V.A.3.-16. - Altitude response during a coordinated tum maneuver by a multi- 
variable McRuer-Krendel pilot inserted in the Black Hawk UH-60A 
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Plot V.A.3.-17. - Pitch angle response during a coordinated tum maneuver by a multi- 
variable McRuer-Krendel pilot inserted in the Black Hawk UH-60A 
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Plot V.B.4.- 1. - Heading response comparison of the OCM pilot and piloted flight data 
during the vertical tracking maneuver. 
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Plot V.B.4.-2. - Rudder pedal deflection response comparison of the OCM pilot and 
piloted flight data during the vertical tracking maneuver. 
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Plot V.B.4.-3. - Altitude response comparison of the OCM pilot and piloted flight data 
during the vertical tracking maneuver. 
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Plot V.B.4.-4. - Vertical rate response comparison of the OCM pilot and piloted flight 
data during the vertical tracking maneuver. 
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TIME (SEC) 

Plot V.B.4.-5. - Throttle deflection response comparison of the OCM pilot and piloted 
flight data during the vertical tracking maneuver. 
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Plot V.B.4.-8. - Longitudinal stick operation comparison of the OCM pilot and piloted 
flight data during the vertical tracking maneuver. 
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Plot V.B.4.-9. - Airspeed response comparison of the OCM pilot and piloted flight data 
during the vertical tracking maneuver. 
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Plot V.B.4.-1 1. - Roll angle response comparison of the OCM pilot and piloted flight 
data during the vertical tracking maneuver. 
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Plot V.B.4.-12. - Lateral stick operation comparison of the OCM pilot and piloted flight 
data during the vertical tracking maneuver. 



Plot V.B.4.-13. - Lateral position response comparison of the OCM pilot and piloted 
flight data during the lateral tracking maneuver. 
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Plot V.B.4.-14. - Lateral velocity response comparison of the OCM pilot and piloted 
flight data during the lateral tracking maneuver. 



Plot V.B.4.-15. - Roll angle response comparison of the OCM pilot and piloted flight 
data during the lateral tracking maneuver. 
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Plot V.B.4.-18. - Rudder pedal operation comparison of the OCM pilot and piloted 
flight data during the lateral tracking maneuver. 
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Plot V.B. 4.-20. - Longitudinal stick operation comparison of the OCM pilot and piloted 
flight data during the lateral tracking maneuver. 
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Plot B. l.-l. - Forward velocity response of the Harrier AV-8B due to a 10 percent 

positive impulse of the throttle. 
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TIME (SEC) 

Plot B.1.-2. - Forward velocity response of the Harrier AV-8B due to a 1 inch impulse 
on the longitudinal stick while in a near hover. 



TIME (SEC) 

Plot B. 1.-3. - Short period response of the Harrier AV-8B forward velocity due to an 

impulse of the lateral stick. 
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Plot B. 1.-4. - Short period response of the Hamer AV-8B forward velocity due to an 

impulse of the rudder pedals. 



Plot B. 1.-5. - Simulation model pitch rate response of the Harrier AV-8B due to a 10 

percent impulse of the throttle. 
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Plot B.1.-10. - Simulation model vertical rate response of the Harrier AV-8B due to a 5 

degree impulse of the nozzle angle. 




Plot B. l.-l 1. - Short period response of the Hairier AV-8B vertical rate due to an 

impulse of the lateral stick. 
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Plot B. 1.-12. - Short period response of the Harrier AV-8B vertical rate due to an 

impulse of the rudder pedals. 
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Plot B. 1.-13. - Pitch rate response of the Harrier AV-8B due to an impulse of the 
longitudinal stick at a near hover. 
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Plot B. 1.-14. - Pitch rate response of the Harrier AV-8B due to an impulse of the 

longitudinal stick at 10 knots. 



Plot B. 1.-15. - Lateral acceleration of the Harrier AV-8B simulation model program 
due to the roll angle of Plot B. 1.-16. 
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Plot B. 1.-16. - Roll angle perturbation of the Harrier AV-8B due to an impulse of the 

lateral stick 
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Plot B.I.- 17. - Lateral acceleration of the Harrier AV-8B due to an impulse of the 

rudder pedals. 
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Plot B.2.-4. - Short period response of the Harrier AV-8B engine speed due to an 

impulse of the throttle. 

2.0 
1.5 
1.0 
0.5 
0.0 
-0.5 

0.0 2.0 4.0 6.0 8.0 10.0 

TIME (SEC) 

Plot B.2.-5. - Short period dynamics of the Harrier AV-8B roll rate due to an impulse 

of the lateral stick at a near hover. 
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Plot B.2.-6. - Short period dynamics of the I’ 

of the lateral stick at 



Plot B.2.-7. - Short period response of the Harrier AV-8B yaw rate due to an impulse 
of the rudder pedals at a near hover. 
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Plot B.2.-8. - Short period response of the Harrier AV-8B pitch rate due to an impulse 

of the rudder pedals at a near hover. 
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